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16195 SELF-FORMED CHANNELS SANDY BEDS 


KEY WORDS: Channel beds; Channel flow; Geometry; Laboratory tests; 
Regime theory; Sands; Secondary flow; Stream flow; Stream meandering 


ABSTRACT: The hydraulic geometry and the local features straight, relatively wide 
self-formed channels composed sandy bed are studied laboratory flumes. 
Introduction one loose and the other rigid bank experiments excludes the 
tendency toward meanders and provides good examples test hydraulic regime 
relations. The experiments support few regime relations derived rationally. The 
transverse bed shapes equilibrium have universal profile regardless initially- 
molded cross section, discharge flow, and sand diameter within the scope the 
experiments. theoretical analysis secondary currents straight channels provides 
reasonable interpretation the local features self-formed channels such the 
longitudinal sand ridges, and this analysis suggests the insignificant effect secondary 
currents the size self-formed channels. 


REFERENCE: Ikeda, Syunsuke, “Self-Formed Straight Channels Sandy Beds,” 


Journal the Hydraulics Division, ASCE, Vol. 107, No. HY4, Proc. Paper 16195, 
April, 1981, pp. 389-406 


16184 KARMAN CONSTANT SEDIMENT-LADEN FLOW 


KEY WORDS: Channels (waterways); Coefficients; Fluid flow; Hydraulics; 
Rivers; Sedimentation; Sediment production; Shear stress; Streams; Velocity 


ABSTRACT: Based the principle energy balance and the effective power 
suspension-load, formula for calculating the Karman Constant sediment-laden 
flow presented. The calculating formula may used practice. According the 
variablity sediment effective power, We, the influence the characteristics 
suspension-load follows: when greater than zero, the value 
decreases with the increase concentration suspension-load (mean value); when 
the value the value increases with the increase (mean value); and 
when the value the variablity has relation with the value 
(mean value) directly. 


REFERENCE: Wang, Shang-yi, “Variation Karman Constant Sediment-Laden 
Flow,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY4, Proc. Paper 
16184, April, 1981, pp. 407-417 


16187 PERFORMACE BAFFLED CHUTE WITH 90° BEND 


KEY WORDS: Aprons; Baffles; Chutes; Chute spillways; Energy 
dissipation; Mathematical models; Scour; Spillways; Splash 


ABSTRACT: The effects 90° bend scour and splash baffled chute have 
been studied 1:12 model. The model design followed recommended United States 
Bureau Reclamation practice except for the square bend. Various combinations 
guide vanes and bend floor elevations were examined 12-ft and 24-ft (3.7-m and 
7.3m) wide chutes carrying cfs/ft (3.7 The simplest, plain, square bend 
caused 2.3 (0.7 more scour depth than does standard, straight chute. The 
additional scour may eliminated using guide vanes lowering the bend floor. 
Sever splash beyond the bend and lower chute walls may require slope protection. 


REFERENCE: Higgins, David T., “Performance Baffled Chute with 90° Bend,” 
Journal the Hydraulics Division, ASCE, Vol. 107, No. HY4, Proc. Paper 16187, 
April, 1981, pp. 419-426 


16191 UNSTABLE TURBULENT FLOW OPEN CHANNELS 


KEY WORDS: Flumes; Friction; Hydraulics; 
Hydrodynamic configurations; Hydrodynamics; Open channel flow; Open 
channels; Supercritical flow; Theories; Turbulent flow; Unsteady flow 


ABSTRACT: Controversy the determination the friction factor, for turbulent 
flow open channels rectangular cross section has existed among previous 
investigators. number them found that the Froude number affects the value the 
friction factor, while others did not find such dependence. this study, the 
theoretical expression for the critical Froude number first derived. Given value 
the Reynolds number, the flow becomes unstable above the critical Froude number. 
This critical Froude number, which appears strongly dependent upon the channel 
width, then introduced equation Rouse calculate the values the 
friction factor under unstable flow conditions. The influence the Froude number 
the friction factor appears important, but the channel width determinating 
factor well. 


REFERENCE: Berlamont, Jean E., and Vanderstappen, Nicole, “Unstable Turbulent 
Flow Open Channels,” Journal the Hydraulics Division, ASCE, Vol. 107, No. 
HY4, Proc. Paper 16191, April, 1981, pp. 427-449 


16206 FLOW THROUGH LARGE SUDDEN EXPANSIONS 


KEY WORDS: Expansion; Flow characteristics; Flow measurement; Flow 
patterns; Reynolds number; Shear flow; Turbulence; Uniform flow 


ABSTRACT: Flow patterns for large, sudden expansions have been obtained 
experimentally. The experimental studies reveal that flow patterns for large expansions 
are highly asymmetric and unsteady. Two distinct types flow have been observed 
which the flow characteristics are essentially mirror images each other. Turbulence 
develops peak values earlier and decays faster with increasing expansion ratios. 
Because the increased turbulent mixing with increased expansion ratio, the mean 
flow pattern becomes steady and uniform more rapidly with increased expansion ratios. 
For the Reynolds number range 0.5 10° 1.0 10°, there Reynolds number 
influence these flow patterns. 


REFERENCE: Mehta, Prabhulal R., “Separated Flow Through Large Sudden 
Expansions,” Journal the Hydraulics Division, ASCE, Vol. 107, No. HY4, Proc. 
Paper 16206, April, 1981, pp. 451-460 


16209 MASS TRANSPORT UNCONFINED AQUIFERS 


KEY WORDS: Aquifers; Dispersion; Experimental data; Finite elements; 
Flow rate; Free surfaces; Ground water; Mass transfer; Mathematical 
models; Water pollution 


ABSTRACT: The weighted residual technique conjunction with the finite element 
method employed formulate numerical models solve flow and transport 
problems nonstationary saturated ground-water zone. Through the formulation 
presented, elements are allowed move accommodate the nonstationary saturated 
domain without the necessity for interpolation when proceeding from one time step 
the next the solution the mass transport equation. The models can simulate the 
flow water and the movement pollutant the saturated zone homogeneous 
isotropic aquifer. The moisture influx from the unsaturated zone can exert some 
influence over mass distribution the neighborhood the free surface. 


REFERENCE: Guvanasen, Varuttamadhira, and Volker, Raymond E., “Simulating 
Mass Transport Unconfined Aquifers,” Journal the Hydraulics Division, ASCE, 
Vol. 107, No. HY4, Proc. Paper 16209, April, 1981, pp. 461-477 
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16213 SHALLOW WATER FLOW OVER POROUS MEDIA 


KEY WORDS: Computation; Drainage; Fluid flow; Ground water; 
Hydrography; Infiltration (water); Mathematical models; 
Numerical analysis; Open channel flow; Open channels; Porous media; 
Runoff; Seepage; Water flow 


ABSTRACT: conjunctive mathematical model for surface-subsurface flow system 
developed. The unsteady surface flow explored set one-dimensional 
dynamic wave equations for shallow water, which are solved using four-point 
implicit finite-difference scheme. The transient subsurface flow two-dimensional, 
with potential gradients vertical direction well the surface flow direction. 
The porous medium may saturated unsaturated both, and can 
nonhomogeneous and anisotropic. The subsurface flow equation solved employing 
successive line over-relaxation implicit finite-difference scheme. The surface and 
subsurface flow components are coupled the ground surface considering the mass 
conservation and pressure relationships. The model verified using existing 
experimental data and analytical solutions for special cases. 


REFERENCE: Akan, Ali Osman, and Yen, Ben Chie, “Mathematical Model 
Shallow Water Flow over Porous Media,” Journal the Hydraulics Division, ASCE, 
Vol. 107, No. HY4, Proc. Paper 16213, April, 1981, pp. 479-494 
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SELF-FORMED STRAIGHT CHANNELS SANDY BEDS 


INTRODUCTION 


Water and sediment alluvial rivers interact yield self-formed channel 
geometry. The simplest but most essential morphological problem self-formed 
channels how interpret rationally the process side-bank erosion and 
the resulting stable channel formation straight channels with noncohesive 
materials. 

The river engineers India established the regime theory the basis 
their experiences obtained excavating and maintaining large canals with straight 
alinement. Any method based dynamics has not succeeded deriving the 
regime theory rationally. For long time, the only analytical model has been 
the threshold theory, which assumes the critical tractive force everywhere 
canal section (see e.g., Refs. and 10). 

Experiments the stable channels coarse noncohesive sands were made 
Wolman and Brush (16), and they showed that the concept static equilibrium 
not realistic. There exists dynamic equilibrium for which the sediment 
actively transported least the central bed region and bed load vanishes 
the side-bank regions. The threshold theory cannot applied most rivers, 
because sediment transport occurs without widening their widths the rivers. 

theory which allows for bed-load transport was first developed Hirano 
(6). treated analytically the process widening straight channels where 
the sand grains are assumed move down the bank the center bed region. 
However, sediment transport seems incompatible with the stable channel 
concept. order for the sediment constituting the bed transported, the 
shear stress must greater than the critical stress the center bed region. 
Since the bank region joins the center bed region smoothly, the shear stress 
must above the critical stress somewhere the bank region, and bank erosion 
will again occur. Thus, long the bed load exists, the channel must widen 
indefinitely. This incompatible with the aforementioned concept graded 

Prof., Dept. Foundation Engrg., Saitama Univ., Urawa, Saitama, Japan. 

Note.—Discussion open until September 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication July 28, 1980. 
This paper part the Journal the Hydraulics Division, Proceedings the American 
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rivers with nocohesive bed materials. This paradox vividly described Parker 
(12), and termed the channel paradox.’’ resolved the paradox 
taking account the lateral momentum transfer due turbulence, and 
succeeded deriving several rational regime relations for straight reaches. The 
relations predict the gross profiles stable rivers, such the depth the 
width slope, discharge, and sediment characteristics are specified. 

These works still leave several important phenomena, such the process 
self-formation channels, the lateral stable bed profile, the effect secondary 
flow the formation channels, and on, explained more fully. 

the research described herein, the process widening and hydraulic 
characteristics self-formed straight channels with noncohesive sands were 
studied two laboratory flumes which are different size. The experiments 
were conducted under idealized condition exclude the effects several 
disturbances such ripples, dunes, and meanders. The secondary circulations 
found the experiments are analyzed theoretically, and are proven 
consistent with the observations. 


Apparatus and Procedure.—Two flumes were used for the experiments. One 
flume (for the sake simplicity termed the long and 
has width cm; the slope and discharge are adjustable. The side walls 
the flume are made glass, which allows the convenience visual 
observations. Water circulated pump, and the discharge measured 


with triangular weir equipped the bottom tank. adjustable tailgate controls 
the flow depth the downstream end, not yield the drawdown backwater. 
well-graded sand with median diameter, 1.3 mm, and the gradation 
coefficient, 1.3 were used, because for this size ripples dunes were 
not produced. The flume was covered the sand the reach between the 
section and 14-m section from the upstream end. The initial cross section 
was trapezoidal and molded template. Only one side bank was made mobile 
exclude the tendency toward meanders. The molded flume was filled with 
still water before the start experiments. Then the tailgate was opened and 
the water was recirculated. 

The elevations free surface and bed were measured point gage 
the 9-m section and section from the entrance where the width and the 
depth formed channels were uniform, and the effects flume ends were 
minimized. Sand feeding was made the entrance prevent degradation 
the upper reach. However, the degradation was not completely avoided, because 
sands were not supplied sufficiently from the side bank the upper reach. 
The sediment transport rate was measured with trap cage every 
along the traverse after stable state was achieved (as revealed later 
figure). the same time the distribution flow velocity cross section 
was determined with pitot tube for several runs. series runs were 
conducted the large flume, and generally required reach stable 
State. 

Small-flume runs were conducted 2-m long, 30-cm wide recirculatory, 
tilting flume. well-graded sand which had median diameter 0.85 mm, 
and the gradation coefficient 1.3, were used. The procedure the experiments 
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was about the same for the large-flume runs. The sediment transport rate 
and the flow velocity were not measured, because the flume was too small. 
Fourteen runs were performed the small flume, and ripples, dunes, 
meanders were found throughout the total set runs. 

View Representative detailed description representative run 
will give general overview the series experiments, because the features 


Elevation in cm 


t in hr 


FIG. 1.—Variation Free Surface Elevation with Time for Run 


c 
N 
~ 
a 


t in hr 


FIG. 2.—Variation Free Surface Width with Time for Run 


the experiments did not differ essentially with each other. Since the total 
process necessary reach stable state was studied detail for Run 17, 
has been selected the representative run. The initially-molded cross section 
was trapezoidal with lateral slope 30°, bottom width, 23.5-cm 
top width, and 7-cm depth. The experiment was continued for total run time 
until dynamically stable state had been achieved. The discharge was 


7.5 
7.0 
6.5 
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4,130 throughout the run. The water surface elevation was initially about 
below the side bank, and rose little with time. The water surface 
height increased about 2.5 during the first and then reached stable 
level shown Fig. The massive fall sands reported Ashida, 
al. (1) did not occur the water margin. The sand particles the margin 
fell one one, and the width increased smoothly with time. stable width 
was achieved after about had passed from the beginning (see Fig. 2). 


Thickness in cm 


t in hr 


FIG. 3.—Variation Sand Layer Thickness with Time for Run and 
y=0 


@ t-Onr 4 te2hr 
Ot-0.5hr @ te4hr 
telhr © t=l2hr 


FIG. 4.—Variation Lateral Bed Topography with Time for Run at: (a) 
(b) Both Sections are Located Lower Reach Large 


Only one particle layer thickness was eroded during the next The eroded 
particles moved down the side-bank slope settle the central bed region, 
and the bed aggraded very rapidly the first hour (see Fig. 3). However, 
the aggradation was not sufficient the reach between the section and 
6-m section despite the sand feeding the upstream end due the shortage 
supply from the side bank. uniform flow depth and uniform bed thickness 


4 
2 
4 
. 0 10 20 30 
y inom 
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were established the middle and lower reach the stream. The sediment 
was actively transported the central bed region and vanished the side-bank 
region after stable state was established. 


TABLE 1.—Summary Experimental Data 


Q/2, 
cubic Mean depth, slope Total 
centimeters B/2, free Reynolds |run time, 
per second centimeters centimeters surface hours 


(5) 


4 


The variation the lateral bed profile with time shown Figs. and 
4(b). The profiles are nearly identical, and concluded that uniform and 
stable reach was established the channel. 

The experimental data are summarized Table 


265 6.20 1.63 1/280 2.5 
268 6.50 1.57 1/200 3,770 
310 7.90 1.46 1/260 3,740 
310 7.25 1.66 1/245 4,200 1.5 
340 7.60 1.68 1/226 3,800 
346 7.05 1.74 1/270 1.5 
350 7.30 1.68 1/310 1.5 
361 7.95 1.64 1/240 4,220 1.5 
408 8.00 1.53 1/310 1.5 
413 7.60 2.07 1/280 4,210 
467 8.70 2.07 1/333 3,200 
530 9.90 1.83 1/490 
530 9.65 2.01 1/280 
630 11.7 2.05 1/240 
3,800 31.0 3.43 1/467 8,400 10.5 
4,070 32.2 3.54 1/479 8,500 
4,130 32.4 3.59 1/466 8,600 
4,480 35.2 3.61 1/506 8,700 
4,550 33.7 3.74 1/550 
4,630 35.7 3.73 1/498 8,700 
4,680 37.5 3.43 1/450 
4,780 35.4 3.79 1/502 8.3 
4,820 37.6 3.39 1/434 8,700 
4,870 35.6 3.76 1/530 11.5 
4,940 37.5 3.61 1/464 9,000 
4,980 39.1 3.55 1/468 8,600 
5,130 39.0 3.61 1/495 8,700 
5,290 39.3 3.70 1/496 9.5 
5,440 39.8 3.78 1/493 9,400 
5,740 39.6 3.95 9,900 
5,990 39.9 3.91 1/508 10,500 8.3 
6,300 41.7 4.02 1/503 10,200 
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Data 


known that the width significantly correlated discharge, and several 
empirical formulas expressed power-law relations have been proposed: 


which width free surface dominant discharge. Variables and 
numerical constants. 

Leopold and Maddock (11) found that the exponent has average value 
0.5 based extensive field data. test whether the relation applicable 
also laboratory flume data, the measured values width are plotted against 
the discharge Fig. Though there exist considerable scatters especially 
flume data Wolman and Brush (16), the aforementioned relation seems 
hold least for each kind sand. The constant however, different 
each kind. The stable width seems less the diameter sand becomes 
large, even discharges flow are identical. The power law therefore not 


Wolman and Brush (0.67mm sand) 
Wolman and Brush (2.0mm sand) 
Present study (0.85mm sand) 
Present study (1.30mm sand) 


c 
N 
~ 
o 


Q/2 in cc/s 


FIG. 5.—Half-Width Plotted Against Half-Discharge Flow 


rational that does not account for any effect bed materials. The stable 
width active sandy channel should determined combination 
several hydraulic variables such the shear stress and the diameter sand 
particles. 

Parker (12,13) obtained rational regime relations for active gravel rivers (ignoring 
secondary flows): 


first regime relation 


0.830 
0.866 (2) second regime relation 


HY4 
10? 
op .* 
os 
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specific gravity sediment and 1.65 for natural sands; 


10? 


Line of perfect agreement 


FIG. 6.—Test Eq. for Laboratory Data; Symbols Fig. 


~ 
c 
> 
> 
2 
o 


Line of perfect agreement 


Predicted B in cm 
FIG. 7.—Test Eq. for Laboratory Data; Symbols Fig. 


The first regime relation tested Fig. with the same data Fig. 
The theoretically predicted values agree fairly well with the laboratory 
data, and the agreement with the large-flume data excellent. 
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The second regime relation examined also for the same data (see Fig. 
7). predicting the stable width, Eq. calculated from the first 
regime relation, because three parameters, and are thought 
independent and primary hydraulic variables alluvial rivers. While considerable 
scatter present, systematic discrepancies between the four kinds sand 
cannot detected. conclusion, the first and second regime relations are 
reasonable, and appear apply also laboratory flumes. 

The laterally averaged transport rate sediment plotted Fig. which 


Run 16 19 26 29 32 
a 0.51 0.38 0.86 1.00 1.24 


relation 


FIG. 8.—Bed-Load Relation for Runs; Written Grams per Minute 
Centimeters 


13 
7 
25 (x=9m) 
25 (x=11m) 
31 


FIG. Side-Bank Profiles for Stable State 


and shear velocity. The measured points correspond low transport 
rates, implying suspended sediment. The third regime relation (bed-load 
relation) was found disagree with the data. 

Local Features Self-Formed Channels.—The stable lateral bed shapes 
the side-bank region (including small-flume runs) observed the present work 
are plotted Fig. which the symbol displacement thickness defined 


which local depth; center depth; and lateral coordinate from 
the water margin. The shapes thus nondimensionalized are roughly similar 


0.0 
Run a 
0.2 
: 
@ Run 
@ Run 0.4 
0.6 
0.8 
1 B/2 
(5) 
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regardless discharge, scale flume, and kind sand. well-expressed 
exponential function defined 


The measured shapes cannot simulated with cosine curve derived 
Glover and Florey (5) under the assumptions static equilibrium and noncohesive 
materials. The surface tension acting between sand particles near the water 
margin seemingly negligible, because any effect cannot detected the 
shapes. 

Fig. provides example dimensionless bed profiles corresponding 
time increments. While considerable scatter exists, such profiles maintain 


t=30min. 
t=60min. 
t=120min. 
t=240min. 
t=650min. 


FIG. 10.—Dimensionless Side-Bank Profiles Corresponding Time Increments for 
Run 


similarity the side-bank region during the widening process, and are also 
well-expressed Eq. The profiles cannot made similar nondimensiona- 


which shows functional relationship. The ratio the central bed region, 
B,, the side-bank region, B,, varies with time the widening process 
(see Fig. 4), and the function varies with time. 

Longitudinal sand ridges were formed several runs (see Fig. 11). With 
crests several millimeters above the bed, the parallel ridges extended unbroken 
from little downstream the entrance the downstream end the movable 
bed. Some the ridges were very stable and persisted almost throughout the 
run time. The distance between the crests was roughly twice the depth, 
reported hitherto (see Refs. and 16). example the bed topography, 
which ridge elevation exaggerated factor two the vertical, 
shown Fig. 11. natural suppose that these ridges may have their 
origin series paired vortices with opposite directions rotation. The 


(-£) 
0.0 
0.8 
1.0 
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scheme the secondary currents inferred from the bed topography shown 
the figure. 

minute observation Fig. reveals that for Runs 13, 17, and 31, where 
the sand ridges were formed, there are systematic discrepancies between the 
exponential profile and the experimental values. The local depths are larger 


FIG. Longitudinal Sand Ridges Formed Run 17; Arrows Indicate 
Ridges which Extend Unbroken From Upstream Downstream End 


y in om 


FIG. Bed Profile Corresponding Fig. 11; Measurements were Made 
and Broken Lines Indicate Secondary Currents Inferred From 
Bed Topography; Note That Vertical Scale Exaggerated 


Center bed regione—O—e Side bank region 
ac 30 40 
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than the solid line around and they are shallower 
where the sand ridges were yielded. For Run 25, which the sand ridges 
were not observed, the experimental values can simulated with sufficient 
accuracy the exponential function. can therefore imagine the existence 
pair vortices also the sloping side-bank region. 

The velocity distributions were measured for Runs 16, 19, 24, 29, and 32. 
The isovels show marked lateral undulations Runs and (see Fig. 13), 
and they are thought due the secondary currents. Indeed, comparison 
the bed topography Run (see Fig. 12) with the isovels Run 
reveals that the longitudinal fluid velocity slow immediately above the ridges 
and swift between them (both runs are identical within the error measurements). 

Corresponding the undulations, the lateral distributions bed-load transport 


show wavy profiles. The bed-load transports are diminished the side-bank 
region for all runs. 


in gr/cm 
° 


FIG. 13.—Lateral Distributions of: (a) Bottom Shear Stress; (b) Sediment Transport 


first estimate the lateral distributions the bottom shear stress can 
obtained from the normal depth method defined (see, e.g., Ref. 15) 


which mass density fluid; area cross-sectional segment 
between adjacent straight lines normal bottom; and wetted perimeter 
corresponding the segment. The calculated bottom shear stress described 
with broken line Fig. 13, which the threshold shear stress for level 
and flat boundaries also estimated with the Shields curve. While the lateral 
inclination the bed naturally affects the incipient motion, the estimation may 
give reasonable values least the central bed region. The bed-load transport 
rates are diminished show abrupt decreases the critical points thus estimated. 
The bottom shear stresses exceed the critical about 20% the major parts 
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the central bed region for the runs, which exactly one the constraints 
imposed applying the regime relations Eqs. and 


described already, the secondary currents seem affect least the local 
features the hydraulic variables such the local depth, the isovels, and 
the sediment transport rate. Analysis the secondary flows should provide 
more rational explanation the local features, especially the sand ridges. 
The analysis will moreover provide useful information for the adequateness 
the rational regime relations (the effect secondary currents neglected). 

Equations Motions.—It can shown that for steady uniform turbulent 
flows straight channels the Reynolds equations motions reduce (see 
Refs. and 


which the pressure terms are eliminated cross differentiation, and the 
convective and molecular diffusion terms are neglected Gessner and Jones 
(4) indicated wind tunnel with square cross section. The symbols and 
the coordinates the transverse and vertical directions, respectively; 
the kinematic Reynolds shear stress; and and the Reynolds normal 
stresses and directions respectively. 

The Reynolds shear stress can equated the temporal mean velocities 
with component the eddy diffusion tensor, 


which and the temporal mean velocity components the and 
directions, respectively. Introduction stream function the secondary 
currents gives simpler form Eq. 


which Using the Brundrett and Baines (2) measurements 
the Reynolds shear stress and the secondary flow, estimated 
crudely identical (see Ref. 


which the Karman constant; and spatially-averaged boundary 
shear stress. the very crude assumption constant eddy viscosity allowed, 
then substitution Eq. into Eq. gives secondary flow equation 


The secondary flows corners square ducts similar pair vortices 
seen the left corner Fig. 11) were examined based Eq. 13, which 
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was shown give reasonable estimation the actual secondary flow field 
(see Refs. and 7). So, can postulated that Eq. interprets the secondary 
flows straight channels least the lowest order. 

Secondary Currents Channels.—As mentioned previously, steady 
uniform secondary currents exist also turbulent flows wide-open channels. 
impossible suppose that the driving force the secondary flow, 
Eq. distributed uniformly along the traverse, and slight fluctuation 
will produce the secondary flow. The distribution along the 
traverse thus can assumed take the form 


which coefficient shown approximately unit Perkins (14); 
amplitude the perturbation, much smaller than unity; dimensionless 
lateral wave number; and dimensionless vertical distribution the 
driving force. Corresponding the aforementioned expression, the stream 
function the secondary currents should take the form 


which dimensionless vertical distribution the stream function, 
determined; and Substituting Eqs. and into Eq. 
yields dimensionless equation 


which z/H. The vertical distribution the driving force can 
approximated (see Ref. 14) 


Eq. can solved under the following constraints: (1) must real; 
the free surface and the bottom. The solution for given 


which and on. apparent that the largest cell corresponding 
predominant, and other smaller cells contribute little The 
higher modes thus can neglected, and the velocity components the secondary 
circulation can approximated 
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Relations with Longitudinal Ridges and Isovels.—The theory thus obtained 
predicts the depth-scale secondary currents with opposite directions rotation. 
Once series such pair vortices set array, sand ridges will formed 
where the flows converge. The sand ridges thus formed will next distort the 
field the Reynolds normal stresses along the traverse, which will, turn, 
maintain the secondary currents. The mutual interaction will thus reach stable 
state. Unfortunately, any functional relationship between the characteristic scales 
the ridges (e.g., height the ridges), and the distortion the Reynolds 


Cos 


cm/s 


c 
= 


FIG. 14.—Lateral Distributions and z/H Flume with Artificial 
Sand Ridges; Measurements were Made Kuroki and Kishi (9) 


stresses not yet established, and further analysis this instability mechanism 
the ridges, then the instability problem would resolved completely. 

The distortion the flow fields with the boundaries simulating the 
sand ridges were studied Kuroki and Kishi (9) tilting flume long 
and 0.8 wide. Spacing artificial ridges interval twice the depth, 
they measured the transverse distribution the longitudinal velocity component, 
the vertical secondary flow component, and the Reynolds shear stress, 
uw. The mean flow fields show relatively regular undulations with wavelength 
the interval the ridges, and the field described fairly well 
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0.2 selected Eq. (see Fig. 14). The velocity component transports 
longitudinal momentum deficit upward the crests, resulting low fluid 
velocity above the crests. The situation exactly the same that found 
Runs and the present study. Moreover, they noticed that the regular 
waviness the field was broken the interval-depth ratio was changed 
from 

Experimental data the distance between the adjacent ridges were presented 
Kinoshita (8) and Wolman and Brush (16). Their data are compared with 
the prediction that the spacing sand ridges should twice the depth (see 
Fig. 15). apparent that the ratio maintains the nearly constant value 
two regardless the Froude number and the Reynolds number (3,300-19, 100). 

Effects General secondary flows are thus cellular and 
composed depth-scale pair vortices the central bed region. The local bed 
topography reveals the existence pair vortices also the side-bank region. 
can therefore conclude qualitatively that the secondary currents exchange 


Theoretical line 


@ Kinoshita 
© Wolman and Brush 


the longitudinal momentum only the scale depth, and the lateral momentum 
transfer insignificant. This qualitative examination supports the assumption 


Parker (12,13), and this may the reason why the regime relations describe 
the actual well. 


The size and the shape self-formed channels uniform noncohesive sands 
were studied two laboratory flumes where the center the channel cross 
sections was replaced vertical side wall; the resulting half-channel helped 
suppress meandering tendencies. 

These experiments revealed the following: (1) Corresponding combination 
the hydraulic variables, stable channel which allows bed-load transport 
the central bed region exists; (2) the transverse bed shapes equilibrium 
have universal one the side-bank region regardless the initial shape, 
the discharge, and the sand diameter. The side-bank region keeps the universal 
shape also the widening process, and well-described exponential 


Fr 
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function, rather than cosine; (3) equilibrium depth and width depend not only 
discharge, but also slope and bed materials; and (4) secondary currents affect 
the local features the self-formed straight channels, but seem not influence 
the general features such the equilibrium depth width. 

crude theoretical argument the secondary currents was made based 
the Reynolds equations motions, and explained fairly well the local 
features the self-formed channels. 
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The following symbols are used this paper: 


oPR wn & 


aspect ratio; 

numerical constant Eq. 

width free surface; 

width side-bank region; 

dimensionless width; 

numerical constant Eq. 

median diameter size; 

segmented area; 

wetted perimeter corresponding dA; 

thickness sand layer center bed region; 
(p,/p) submerged specific gravity sand; 
functional relationship Eq. 

dimensionless vertical distribution ww; 
gravitational acceleration; 

center depth; 

local depth; 

dimensionless vertical wave number cell; 
positive integer; 

water discharge; 

volumetric sediment discharge; 

dimensionless water discharge; 

dimensionless sediment discharge; 

sediment discharge per unit width; 

dimensionless sediment discharge per unit width; 
dimensionless depth; 

water surface slope; 

time; 

longitudinal, lateral, and vertical temporal mean velocity, respec- 
tively; 

laterally averaged shear velocity; 

the kinematic Reynolds stress; 

the kinematic Reynolds normal stress; 
longitudinal coordinate from the entrance flume; 
lateral coordinate from the center channel; 
vertical coordinate; 

z/H dimensionless vertical coordinate; 
numerical constant Eq. 14; 

displacement thickness side bank profile; 
dimensionless amplitude perturbation Eq. 14; 
component eddy diffusion tensor; 

the Karman constant; 

lateral coordinate from water margin; 
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water density; 

sediment density; 

bottom shear stress; 

critical bottom shear stress (Shields); 

dimensionless bottom shear stress; 

stream function secondary flow; 

dimensionless stream function; and 

vertical distribution dimensionless stream function. 


406 HY4 

p => 

= 


APRIL 1981 


JOURNAL THE 
HYDRAULICS DIVISION 


VARIATION KARMAN CONSTANT 
SEDIMENT-LADEN FLOW 


Shang-yi 


INTRODUCTION 


The study the law velocity distribution open-channel flow carrying 
sediment load one the basic problems the exploration river bed 
development. Practice has shown that velocity distribution along the depth 
sediment-laden open-channel supposed follow approximately the logarithmic 
formula velocity distribution. the formula there very important factor 
(K) called the von Karman coefficient. vast amount research works have 
been done different workers the study the value 
11,12,14). 

Besides the characteristics suspended load, the factors which have 
influence the variation the coefficient should also include the bottom 
regime stream flow. the study the influence suspension-load 
characteristics the value many researchers had attempted conduct 
experiments under controlled bottom conditions. Owing the complexity 
the problem, however, such attempts have not been successful. 

this paper, based the energy balance movable bed model without 
suspension load and the principle effective power suspension load, the 
following formula presented for calculating the Karman constant the 
region main flow open channel carrying sediment load given 

Where the first term the right-hand side the equation expresses the 
influence the bottom condition the value and the second term 
expresses the influence characteristics suspension load Through 
analysis data the laboratory and natural open-channel flow, value 
the coefficient together with method for calculating the numerical value 
K,, has finally been proposed. Equation calculates the value which 
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Note.—Discussion open until September 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication June 1980. 
This paper part the Journal the Hydraulics Division, Proceedings the American 


Society Civil Engineers, Vol. 107, No. HY4, April, 1981. ISSN 0044- 


407 


408 APRIL 1981 HY4 


given, and verified reasonably well with some data observed both the 
laboratory and the field. 


ANALYSIS 


From the viewpoint energy balance (1), may shown that the energy 
expended support the sediment load suspension not taken directly from 
the effective potential energy, but only part the total frictional loss 
stream flow. 

The velocity distribution reflects intrinsically the resistance drag stream 
flow. For open channels, the resistance drag not only depends the characteris- 
tics sediment load, but also more closely the bottom which 
consist the bed form and the characteristics the bed load (2). must 
consider, therefore, the influence bottom conditions the stream flow when 
studying the effect the suspension load the velocity distribution. Owing 
the complexity bed material motion which acted the bottom 
regime directly, extremely difficult rule out the influence bottom 
conditions throughout the experiment even the bed form fixed and the 
sediment controlled without siltation. 

order facilitate analysis, considering the sediment-laden open-channel 
flow whole model suggested which the effect suspension load 
the velocity distribution the region main flow differentiated from 
that caused bottom conditions. 

Suppose that model open-channel flow with suspension load but 


with the same movable bottom conditions there depth flow and hydraulic 
gradient actual sediment-laden open channel. may then write: 


which the velocity arbitrary point measured from the water 
surface sediment-laden open-channel flow; and and U,, the shearing 
stress and velocity the corresponding point the model, respectively. 

Based the principle energy balance sediment-laden open-channel flow, 
caused the difference velocities due the suspension load and induced 
the effective potential energy supplied the stream flow, while the other 
part caused the differential shearing stresses. The term 
the corresponding differential, part which caused the difference 
velocities and part the local energy dissipated while the other part 
the energy dissipated caused the difference shearing stresses; the 
term U,,)] the corresponding differential, part which caused 
the difference velocities and induced from the transmitted energy, 
and for the same reason, the other part transmitted energy also caused 
the difference shearing stresses. the previous analysis, all the higher 
order differential parts are being neglected. 

the sense the aforementioned energy stream flow and the principle 
effective power sediment, three suppositions could mentioned follows: 

First Supposition.—The changed portion the effective potential energy caused 
the difference velocities arbitrary point proportional the effective 
potential energy supplied the suspension load. may put 
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which proportional coefficient; the term (y, y,,) CUJ denotes 
the effective potential energy supplied suspensionload per unit volume 
arbitrary point per second; the concentration specific volume 
the arbitrary point; the velocity the arbitrary point; the hydraulic 
gradient; and and the specific weights sediment and clean water, 
respectively. 

Second Supposition.—The transmitted energy previously 
mentioned proportional the power suspension load. may put 


suspension load per unit volume arbitrary point; and the settling 
velocity sediment still water. 


Substituting Eq. and Eq. into Eq. obtain 


Putting y,, and substituting into Eq. obtain 


Integrating the aforementioned equation along the depth turbulent region 
stream flow, may have 


the depth stream flow; and the thickness bed layer. 

Third Supposition.—The proportional coefficient equal and they 
have constant value along the depth stream flow. may put 


Substituting into Eq. may obtain 


J 


which the term UJ) denoted Eq. the effective power sediment 
(17); and the depth main stream flow. 


Ordinarily, may write the equation velocity distribution along the depth 
stream flow the following: 


which and the von Karman constant two kinds channels 


(2) 
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previously mentioned, respectively; U,, the frictional velocity; and 
the mean velocities along the depth stream flow the two kinds 
channels previously mentioned, respectively: 


H-Z 


where correction factor; and the mean values the corresponding 
factors along the depth stream flow; and denotes the mean value 
the effective power sediment along the depth stream flow. Setting 


and substituting Eqs. and into Eq. 13, then 
Vw U,, J 


denotes the variation von Karman constant the region 
main stream flow open channels. According the analysis previously 
mentioned, the term 1/K,, part the influence caused the bottom 
conditions the velocity distribution. The second term which lies the right-hand 
side Eq. caused the characteristics suspended load. 

According the variability effective power sediment and from Eq. 
17, three conclusions can drawn follows: (1) When the 
the variability has not relation with the value directly. 


the calculating the value using Eq. 17, necessary that the 
coefficient and determined beforehand. 

three variables a,, o,, and A,. almost impossible that the value 
can obtained theoretical analysis. When want obtain the value 


2.3 U,, 
H-Z 
(12) 
Substituting into Eq. have 
H-h 
(14) 
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Eq. directly, necessary also know the value K,, which 
yet unknown. practically impossible make model open channel 
having exactly the same values and the bottom conditions with the 
prototype, approximate methods must employed solve the problem. 

Ref. may found that there are some experiment data which show 
the characteristics the two kinds sediment-laden open channels which 
the difference between their depths stream flow (AH 0.1 cm) between 
their hydraulic gradients (AJ 0.2 %o), the formation and the form bed, 
and the physical characteristics sediment are approximately the same. 
our analysis, the values are considered the same the two sediment-laden 


TABLE 1.—Calculating Results Coefficient 


centi- 
meters 


(2) 


asa 
per- 
centage 


(4) 


Bed 
pattern 
(8) 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 
Smooth 


Yw U,,J 
The subscripts and denote the corresponding factors two similar 
channels previously mentioned, respectively, and the mean value hydraulic 
gradient equal J,). Using Eq. and calculating with the 
data Ref. may obtain the results shown Table From these 


results, may seen that the values are invariable and may take 
its mean value 


centi- 
per 
sec- per 
(1) (3) (5) (6) (7) (9) 
7.84 2.06 0.261 0.223 0.946 73.5 
7.77 2.06 0.305 0.219 72.8 1.03 
7.74 2.59 0.0192 0.345 1.98 81.8 
7.70 2.57 0.0532 0.330 1.98 79.9 1.41 
7.74 2.59 0.0192 0.345 1.98 81.1 
7.70 2.57 0.0920 0.319 1.98 79.9 1.21 
7.74 2.59 0.0192 0.345 1.98 81.1 
7.70 2.59 0.1230 0.310 1.98 80.5 1.16 
7.70 2.57 0.0532 0.330 1.98 79.9 
7.70 2.57 0.092 0.319 1.98 79.9 1.00 
7.70 2.57 0.0532 0.330 1.98 79.9 
7.70 2.59 0.1230 0.310 1.98 80.5 1.05 
7.70 2.59 0.092 0.319 1.98 79.9 
7.70 2.59 0.1230 0.310 1.98 80.5 1.02 
open channels previously mentioned; may obtain 
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According Eq. 20, have calculated with set practical data which 
has runs and vary the following regions, such the depths stream 
flow, 7.0 195 cm; velocities, 250 cm/s; hydraulic gradient, 


Value 


Value 


(1.2 255) %o; the concentrations sediment, 0.02 0.396 and 
the mid-diameters sediment, 0.018 1.3 mm. 


calculate the settling velocity, use the following equation collective 
system (16): 


which the mean settling velocity without concerning the influence 
concentration sediment along the depth stream flow. The value 
determined the following equation (17): 
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Determination ,,.—Substituting Eq. into Eq. 17, then obtain 
0.26 
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0.755 0.222 log 
calculating and analyzing the practical data, the following results may obtained. 


When 17, the relationship between the value and the 
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FIG. 3.—Regression Results Equation 
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Measured value 


0.44 0.52 
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FIG. 4.—Calculated Results for Data References and 
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(Fig. 1). When 17, K,, decreases increases. the conditions 
U/U,, 24, K,, varies from 0.4-0.2 widely, and until U/U,, 24, 
approaches constant value 0.2 (Fig. 2). 


For the purpose contract, the value clean-water stream flow with 
fixed bed shown also Fig. 


Calculating Equation K.—According the relationship Eq. 20, may 


K K,, Yw 


This the equation for calculating the value and its standard deviation 
about +0.025 (Fig. 3). 


Equation 


have used the experimental data Refs. and for verifying Eq. 
23. The acquired results are shown Fig. and Fig. 


Summary 


Besides the characteristics suspended load, the factors which have 
influence the variability Karman constant may the bottom regime 
stream flow. Ordinarily, the later more important than the former, 
pointed out Vanani and Nomicos Ref. 

The first term the right-hand side Eq. expresses the influence 
the bottom conditions the value and the second term expresses 
the influence characteristics suspended load 

According the variability the effective power sediment may 
conclude the influence the characteristics suspension load 
the following: (1) When the value decreases with the 
with the value directly. 

The calculation for Eq. may used practice. Its standard deviation 
about +0.025. The primary results calculate with the data Refs. 
and are shown Fig. and Fig. 
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The following symbols are used this paper: 


1/(0.755 0.222 log 

specific volume arbitrary point and its mean 
value along depth stream flow, respectively; 
mid-diameter sediment, millimeters; 
depth turbulent region stream flow; 
depth stream flow; 
hydraulic gradient its mean value two, and J,, respec- 
tively; 
Karman constant; 
part determined bottom regime stream flow; 
velocity arbitrary point and its mean value along depth 
stream flow, respectively; 
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velocity arbitrary point and its mean value along depth 
modeling stream flow, respectively; 

frictional velocity; 

settling velocity sediment still water arbitrary point 
and its mean value along depth stream flow, respectively; 
three proportional coefficients; 

specific weight sediment and clear water, respectively; 
thickness bed layer; 

correction factor; 

standard deviation; and 
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PERFORMANCE BAFFLED CHUTE WITH 90° BEND 


INTRODUCTION 


The baffled chute spillway-type structure that carries water lower 
elevations with relatively low terminal velocity and need for constructed 
stilling basin. The Water and Power Resources Service (formerly the United 
States Bureau Reclamation and referred this paper the USBR) operates 
many such chutes its irrigation projects. has developed general guidelines 
for their design provided they are straight plan. However not difficult 
imagine situations where topography interfering structures would make 
baffled chute with one more horizontal bends desirable. 

Questions about the performance chute with bend led the exploratory 
studies reported this paper. The results, mostly about scour and splash, may 
provide design engineer with sufficient information decide whether such 
structure practical solution his problem and whether further model 
tests are order. 


Numerical information about the model will stated terms prototype 
magnitude wherever reasonable so. The model shown Fig. 
included reach approach channel, upper chute, bend, lower 
chute and reach outlet channel. Details the upper and lower chutes 
followed the recommendations the USBR (2) for 24-ft (7.3-m) wide chute 
carrying cfs/ft (3.7 scale 1:12. The discharge and scale 
were similar those straight chute models reported the USBR (1). 

Outside the bend, the model geometry was held constant except for changing 
the width (3.7 four tests and dropping the approach channel 

Assoc. Prof., Dept. Civ. and Environmental Engrg., Washington State Univ., Pullman, 
Wash. 99164. 
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floor (1.5 two tests. Discharge for all tests was (3.7 
Baffles were high and were arranged five rows the upper 
chute and seven rows the longer, lower chute. sill with crest elevation 
zero controlled tailwater stage and sand bed elevation. was located 
(22 downstream the lower chute exit. 

Model Calibration.—The Froude law was used determine model discharge. 
But scaling scour different matter. For these exploratory tests seemed 
appropriate relate bend-model scour the scour tests the USBR (1). 


FIG. Pattern below Straight Baffled Chute (Test No. 0.3 


The model was temporarily modified operate straight chute with 
rows baffles. Several locally available sands were tested order find 
one that would scour about the same way that reported the USBR. 
Results the tests are shown Fig. The scour behavior the river sand 
was sufficiently close that the USBR proceed with the bend-model 
tests. The anomaly between scour depths and relative probably due 
the inability the simple correlate uniquely with complex scour phenomena. 


— 
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Test 


Scope.—The test program was devoted mostly examining scour depth 
function geometry within the bend. Guide vane arrangement and bend 
floor elevation were the principle variables. few tests looked surge height, 
splash against and over the chute walls, and pressures the bend walls. 

each scour test the sand lost during the previous test 
was replaced and the entire sand bed surface, including that within the lower 
end the chute, was brought elevation zero. The model was then operated 


FIG. 4.—Scour Pattern below Baffled Chute with 90° Bend (Test No. 


for min (model time), drained, and the scoured area marked with contour 
lines 1-ft (0.3-m) intervals. 

Some interesting features the scour behavior were noted during these tests. 
Scour patterns closely duplicated themselves under identical test conditions. 
test duration increased beyond min the region bed degradation extended 
farther downstream, but the scour depth the chute exit remained constant. 
the scour hole was deepened hand after 30-min run, few more minutes’ 
running time saw the scour hole refilled its 30-min equilibrium level. 
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More water splashes over the walls the bend model than straight chute 
models. record how high the splash went, chute and bend walls were extended 
vertically adding panels brown chipboard against which splash marks 
showed vividly. record lateral splash extent, chute wall height was held 
(2.7 normal the chute floor and chipboard panels were installed 
that height outboard the chute and parallel the chute floor. 

the event that guide vanes might practical design feature chute 
bends, orifices for pressure transducers were installed the wall the 12-ft 


PLAN VIEW ELEV. OF MAXIMUM 
OF BEND APPROACH SCOUR 
12 24FT REMARKS 
(FT) 


(6) 
12-FT WIDTH 
24-FT WIDTH 


SIMPLE BEND 


t DIRECTION 


w 
ww 


THIN SHEET 
METAL VANES 
11-FT HIGH 


0.75-FT THICK 


1.5-FT THICK 
GUIDE VANE 
8-FT HIGH 


STANDARD PIERS 


7.0 

7.5 

7.0 

7.6 
| 


N 


(3.7-m) wide chute bend. The 12-ft (3.7-m) bend wall would experience about 
the same surging pressure environment guide vane wider bend. 


Tests 


Scour.—The effect the bend scour make deeper and asymmetric 
with respect the lower chute centerline. This effect shown Figs. 
and which, respectively, show scour contours below straight chute (test 
and chute with bend (test 6). 


NO. 
2 
8 IN 12-FT 
CHUTE 
8 
6 8 9.0 
7 6 8.5 SIMPLE BEND 
IN 24-FT 
8 4 8.0 CHUTE 
9 2 7.0 
13 
GUIDE VANES 
8-FT HIGH 
— 
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The complete set scour test results summarized Fig. The additional 
scour depth due the bend varied from 0.3 (0.1 m)-3.6 (1.1 and 
was influenced largely variations geometry within the bend itself. The 
effect entrance velocity the upper chute shown minor tests 
and 

Tests and (Fig. show that square, 90° bend without guide 
vanes 12-ft (3.7-m) chute causes less than 0.5 (0.2 additional scour. 
Tests 6-16 the 24-ft (7.3-m) chute show additional scour depths from 
0.3 (0.1 m)-3.6 (1.1 m). 

Tests 6-9 demonstrate the effect varying floor elevation plain, square 
bend. The lower the floor, the less the scour. But for some reason this trend 
does not exist when guide vanes are present. See tests and 11, and and 
13. 

The curved sheet metal vanes tests and were not practical prototype 
alternative. They were installed simply limiting case flow guidance. 


Intersection bend Envelope 
walls at outer corner completely wet area 


w 


Envelope surge 


Pressure 
orifices 
Be 


> 
2 
WwW 


FIG. and Splash Profiles along Outer Walls 12-ft Chute and Bend 
0.3 


The baffles test were thought practical alternative and, surprisingly, 
caused the deepest scour all. 

the several arrangements plane, segmented vanes (tests 12-15) the single 
vane test may the most practical. 

However, all consideration vane patterns largely academic. Their con- 
tribution scour reduction would probably not repay their cost. Instead the 
simple bend designs tests and would appropriate design choices. 
additional row baffles the bottom the chute would protect against 
the added scour. 

Splash.—The agitated flow the upper chute around and over the baffles 
prolific source splashing water. The trajectories most the drops 
are parallel the chute center line. the usual straight chute these drops 
will land harmlessly farther down the chute the stream below. But 
the bend model significant quantity water thrown beyond the bend 


= 
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shown Figs. and prototype splash all like the splash the 
model, then splash and its overbank erosion potential become item concern 
for the designer. 

Surge.—On tests and the usual stage the choppy water surface 
the bend was about elevation (4.9 m). Surges along the wall rose about 
(0.7 above this level. Fig. shows the envelope maximum surge 
height along the bend and lower chute walls for the 12-ft (3.7-m) chute. Surges 
were bit lower the 24-ft (7.3-m) chute. The surge height envelope along 
the right wall the lower chute was about 0.5 (0.2 lower than along 
the outer wall. 

Pressure.—Pressures against the 12-ft (3.7-m) chute bend wall may good 
approximation pressures acting guide vanes wider bends. Pressures 
were recorded the five locations shown Fig. for the bend test 
They suggest that the lower (2.9 vane will subjected uniform 


Completely wet area 


FIG. 7.—Lateral Distribution Completely Wet Area Left Chute 
0.3 


unbalanced thrust (1.2 water. This includes the effect flow 
separation behind the vane. The frequency the oscillating pressure against 
the bend wall the model was about Hz. 


Summary 


The straight, baffled chute has proven useful hydraulic structure. 
The incorporation horizontal bends would widen its field application 
performance were not seriously degraded. This empirical study explored the 
consequences single, 90° bend chutes and widths (3.7 
and 7.3 discharging cfs/ft (3.7 

Except for details flow guidance within the bend, the model design followed 
USBR recommendations for straight, baffled chutes. The effects several 


Slope 
irection 
= == = information on splash 
[slo igh of chute) 
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guidance schemes scour were measured. few splash tests were also made. 
The results may summarized follows: 


simple, square, 90° bend standard baffled chute carrying cfs/ft 
(3.7 increases the scour depth 0.5 (0.2 12-ft (3.7-m) wide 
chute, and 2.3 (0.7 24-ft (7.3-m) wide chute. 

The additional scour below the 24-ft (7.3-m) chute may virtually eliminated 
placing guide vanes the bend lowering the bend floor. Alternatively, 
one more row baffles may installed the bottom the lower chute 
protect against the additional scour. 

Splash beyond the bend and lower chute walls may require slope protection. 


From these tests the writer concludes that for flow rates least 
cfs/ft (3.7 there important reason based performance why 
designer should not consider baffled chute with bend where could 
solve his problem. 
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UNSTABLE TURBULENT FLOW OPEN CHANNELS 


INTRODUCTION 


review the literature dealing with the determination the friction factor 
for turbulent flow rectangular open channels, revealed controversial fact; 
number authors found that the Froude number affects the value the 
friction factor, while others asserted that may not important parameter 
after all. 

this study, the critical Froude number first theoretically derived. Given 
value the Reynolds number, the flow for the Froude number above the 
critical value becomes unstable with the appearance roll waves. The critical 
Froude number tends strongly dependent the channel width. Then, 
the value the friction factor under unstable turbulent flow conditions can 
calculated using empirical formula developed Rouse. Formulas 
for the friction factor function the Froude number, the Reynolds number, 
and the channel width can thus developed for smooth channels. The present 
study will clarify the apparently contradictory findings different experimenters. 


Equations Motion 
The general two-dimensional Reynolds equations for turbulent flow open 
channels, taking into account the vertical velocity and acceleration, are: 


yy 
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which the x-axis chosen along the channel bottom and the y-axis 
velocity component the direction; slope angle the channel; 
mass density the fluid; and Reynolds stresses. 


Continuity Equation 


Integrating Eq. over the cross-sectional area, with the help the kinematic 
surface condition gives the continuity equation for the mean flow 


which U(x,t) mean velocity over the cross section time and 


Eqs. and may integrated over cross section assuming the following. 


The flow gradually varied, i.e., /ax 
The expression for the mean shear stress the channel walls is, even 


unsteady flow: 
which the Darcy-Weisbach friction factor. 
The velocity profiles are similar all cross sections: 
u(x,y,t) F(z) 


which y/H. This assumption valid both the case laminar flow 
(parabolic profile) and turbulent flow (logarithmic profile). From 

Products lower derivatives and with respect and are 
small, compared with the higher derivatives (11). 


From and follows that 
Eq. becomes, with assumption 


Substituting Eq. into Eq. with the help Eq. and assumption one 
obtains, after integrating Eq. between and 
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P Ty 2 H z 


aH H z * oF 


After integrating Eq. over the cross section, substituting Eq. and taking 
assumption into account, one obtains the momentum equation: 


+o 


which the hydraulic radius; and C,, C,, and constants, the value 
which depends upon the shape the velocity profile F(z), and defined 


z 


Fdz+ 


0 


z * oF 
0 


Dropping the third derivatives, Eq. reduces the equation. 
then reduces the hydrostatic law. 

the contrary, assuming uniform velocity distribution [F(z) 
and 2], but taking the higher derivatives into account, 
one finds the momentum equation, obtained Keulegan (11). Eq. then reduces 
to: 


The coefficients C,, C,, and can calculated for different values 
the turbulent velocity profile approximated one: 
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F(z)= 


Since F(z)dz=1, one obtains: 
0 

1 


Furthermore 


—— 


Parabolic 
Uniform 
Trapezoidal 


Triangular 


Typical values the constants C,, C,, and for various velocity profile 
assumptions are listed Table 


When searching for stability criterion, one should find out under which 
conditions infinitesimal small amplitude disturbance wave will grow decay 
Thus 

au’ (U,+ U’) 


HY4 
(16) 

2-5 
2b) 
(2-5) 
TABLE 1.—Typical Values Constants C,, C,, and 

(1) (4) (5) (6) 

1.20 0.825 2.496 1.864 

0.088 1.03 0.977 2.068 1.094 

0.145 1.05 0.961 2.113 1.161 

0.281 1.10 0.921 2.227 1.341 

1.33 0.750 2.700 2.400 


HY4 UNSTABLE TURBULENT FLOW 
aH’ 
Neglecting the product H’, Eq. becomes: 


aH’ 


Differentiating Eq. with respect and neglecting products and powers 
and greater than one obtains, with the help assumption and 
23: 


+ 2 
pR, 


Applying Taylor series and R(H), and dropping the higher order 
terms, the right side the equal sign Eq. may written as: 


introducing Eq. 23, one obtains: 
(=) aH’ 
(=) aH’ 
pR,H, 


The linear differential equation, Eq. 24, has solution the form: 
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c,; and the infinitesimal small amplitude the disturbance 
wave. 


Introducing Eq. Eq. 24, one obtains, when putting both the real and 
the imaginary part equal zero: 


au/, 


Introducing the dimensionless variables: 


and dimensionless wave celerity; 


the dimensionless wave number; 


the dimensionless wave length, 
0 


3 


(aR 


and 


HY4 UNSTABLE TURBULENT FLOW 
Since, with assumption 


8 au 8 au Pp 4 0» 


The friction losses are best taken into account using the White-Colebrook- 
Thijsse formula (18): 


which k,/R the equivalent relative roughness the channel; the 
hydraulic radius; R/v), the Reynolds number; and the kinematic 
viscosity the liquid. 


TABLE and for Different Friction Formulas 


Friction formula 
(1) 

White-Colebrook-Thijsse 

(Eq. 31) 
Laminar: 64/R 


(3) 


—fo/Ro (@R/dH), 


smooth channel (k, 0), Eq. reduces to: 


The expressions for and are summarized Table for 
different friction formulas. 


Using the White-Colebrook-Thijsse formula, function introduced: 


—— + 
For hydraulic smooth surfaces: 


The stability criterion obtained putting Eqs. and 29. 
yields the critical Froude number, and the celerity, infinitesimal 
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small amplitude disturbance wave F.. The flow will unstable 


0). Then the amplitude the disturbance wave increase with 
time. 


Taking Eq. into account, one obtains: 


(. = —2a Ci. + Cre (Cs C,C,.+ 


c 


which the wetted perimeter. 


TABLE 3.—Values Critical Dimensionless Wave Celerity, for Different Friction 
Formulas 


Friction formula 
(1) 


Chézy (M/2) 
Laminar 1+2M 


Since for rectangular section width 


reduces to: 


For channel infinite width, The value depends upon the 
friction formula chosen and computed from Eq. given Table 

The critical Froude number, F., can obtained from Eq. 35, substituting 
the results Table One obtains, using the White-Colebrook formula: 


(2) 


HY4 


seen that the critical Froude number depends upon: (1) The shape 
the velocity profile (a, C,, C,, (2) the Reynolds number, and parameter 


formula 
a=1 - Manning formula 
formula 
formula 


0 1 


(B=0) (B=a@) 


FIG. 1.—Critical Froude Number, Function for Different Values 
According Resistance Laws Chézy and Manning 


signifying the relative roughness, (3) parameter representing the channel 
width, and (4) the dimensionless wave number, For disturbance 
with infinite wave length), one finds: 


Using the Chézy formula 
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With the Manning formula 


2 2 2 1/2 


(10). 


formula. seen that instability more likely occur wide section. 


FIG. 2.—Critical Froude Number, F_, Function Reynolds Number, and Channel 
Width, for Smooth Channel and 


Using the White-Colebrook formula, Eq. 39, for becomes 
= 


Since its derivation was general, Eq. indeed holds true for /aminar flows 
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which 1.2 (see Table for parabolic velocity profile). One obtains: 


found Ishihara (4). 

Yih (22), using potential theory and the Orr-Sommerfeld equations found 
(not taking into account the capillarity effect): 0.527 and and 
3.33 cotg Assuming logarithmic velocity profile, Iwasa has shown 
(6) that: 


ec 


Fig. shows the critical Froude number, function and the channel 


| Experimental results 


| @ unstable flow 


{ stable flow 
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FIG. 3.—Experimental and Computed Values Reynolds and Froude Number for 


width, computed with Eq. for smooth channel and and both 
for and computed according Eq. Iwasa. For given value 
one obtains marked increase for narrow channels, compared 
with channel infinite width. should noted that, when scale model 
built constant stable flow the model does not necessarily imply 
the same for the prototype. 

Fig. some experimental results are shown for 0.24 (10 ft) and 
smooth channel (1). For this particular value curves constant can 
drawn the (R,F) plane. seen from Fig. that for given channel 
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width, minimum slope exists below which flow instability occurs. For 
given discharge, the flow turns from stable unstable regime 
critical Froude number obtained increasing the channel slope. For fixed 
channel slope, instability may disappear for large values the discharge 
(or R). 

The experimental results not match the theoretical results completely. 
This can explained (see Fig. the fact that slight increase channel 
roughness may increase considerably (although the channel was reported 
smooth, might have been the order few tenths millimeters), 
but especially because takes some distance along the channel for developing 
roll waves (12,1), especially for the greater discharges and steeper slopes, 
that flow conditions reported being stable, i.e., with (visual) development 


channel roughness ks/R 
1. = 0.0000 

2 = 0.0001 

3 = 0.0005 
= 00010 
= 00050 
= 00100 


R 5.105 


FIG. 4.—Influence Channel Roughness Critical Froude Number, According 


roll waves, might yet have shown roll waves the channel had been longer 
than (40 ft) (1). 


Friction Factor 


According dimensional analysis Rouse (15), the friction factor, 
should function not only the Reynolds number, the relative channel 
roughness, and the shape the flow section, but also the Froude number. 
Since the few indications gravitational influence found the literature 
were all obtained Froude numbers well excess unity, Rouse assumed 
that the phenomenon was some way connected with the surface instability 
involved the formation roll waves. 

From experimental results obtained (15,16), and experimental evidence 


| 
-4+ 4+ + +-4 + 


UNSTABLE TURBULENT FLOW 
obtained Nemec (2), Rouse found for smooth channels: 


2.03 log for F>F. 


except for small difference the constant value, which agrees with the laws 
Blasius, Thijsse (18), Keulegan (9), Straub (17), and others. 


FIG. Function and Froude number, for Different Values 
Channel Width, 


Chen (2), Powell (13), and others wrote 


which and are some function and for and constants 
for 

Introducing the value (see Eq. 40) into Eq. Rouse, one finds 
that should not only function and but also marked function 
the channel width. The value computed with and 42, 
All stable flow situations are situated straight line the 
plane (see Eq. 32). Similar graphs can drawn with sin 
parameter instead From Fig. can seen that the influence 

Since for smooth channel, and given channel width, function 
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only (see Eq. 40), Eq. can approximately written 


linear regression computed results similar those presented Fig. 
yields, for the best straight line, according Eq. 44, the values and 
given Table valid for 0.014 0.040. From Table can seen 
that nearly independent the channel 


TABLE 4.—Results Linear Regression Analysis 


Correlation factor, 
percentage 


(5) 


Eq. 42. 
Eq. 32. 


obtain the value the friction factor hydraulically rough channels 
Froude numbers excess F., Eqs. and should used together 
with the formula Rouse (16) 


= 
Vf ane 


which the roughness concentration. 


Comparison 


Experimental verification the influence the Froude number the value 
the friction factor, hydraulically smooth channels, was obtained 
Nemec (2,15), Jegorov (8), Powell (13), Iwagaki (5), and Wakhlu (21). Other 
experimenters (19), however, did not find any influence the Froude number, 


7 

(1) 

0.08 2.32 99.75 
(0.3) 

0.16 2.24 99.51 
(0.5) 

0.24 2.21 99.51 

0.40 2.17 99.64 
(1.3) 

1.00 2.12 99.80 
(3.3) 
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attributed the measured differences errors the measurements, both 
(2) (see Table 5). 


TABLE 5.—Experimental Results Found Literature 


Influence 


Reynolds Froude Froude 
Reference meters (feet) interval interval 


(1) (3) (4) 
(8) 10° 5.6 0.9 3.5 


(13,14) 10* 10° 2.49 3.28 
(5) 
5.2 10* 6.7 10° 
2.5 1.5 10° 0.50 12.5 


8.8 10* 2.4 10° 1.19 2.88 


Experimental results by Jegorov 
(B= 0.031 

F=15 

@re25 


—— Authors’ theoretical 


5.107 av 


FIG. Results Obtained Jegorov Comparison with Computed 
Values 


Since the width the various experimental channels (2,5,8,13,14,15,19,21) 
varied between 0.03 (1.2 in.) and 1.10 (3.6 ft), may true that the 
difference channel width might partly account for the contradictory results 
the influence the Froude number the friction factor. 

Jegorov (8) found increase 10%-20% for 0.9 3.5. For lower 
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values the obtained values were identical those obtained from 
the Blasius formula. For 3.5, decreases the Blasius value. The maximum 
increase was obtained for 2-2.3. Some measured values Jegorov 
1.5 and 2.5 are compared with our computed values. All values 
2.5 are close the curves and (see Fig. 6). The influence 
was explained Jegorov (8) studying the stability the interface layer 
between the laminar boundary layer and the turbulent flow. Powell found (14) 
experimentally: 


Eq. according Brock, and reworking Powell’s experimental 
data (2). 


lwagaki 


—— Authors’ theoretical 
(B=0.40m) 


10° 


FIG. 7.—Results Obtained Comparison with Computed Values 


Eq. was obtained from experiments with (see Table 5). increase 
decreasing value with increasing which variance with Eq. 44. 
may that the narrow interval values did not allow the extrapolation 

Iwagaki (5) found 0.89 and 
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the last term being small, compared with the others. shown Fig. 
Eq. yields larger values than Eq. 42. 

Wakhlu (21) also found increase with He, too, obtained larger 
values than the theoretical ones obtained means Eq. 42. 


Tracy and Lester's experimental 


© 100 < F< 125 
#125 «175 
@175 225 
© 225 < 275 
@ 325<F 


—- Tracy and Lester's theoretical 
— Authors’ theoretical 


(B=107 m) 


5.103 10° 


FIG. 8.—Results Obtained Tracy and Lester Comparison with Computed Values 


Tracy and Lester (19) found for both subcritical and supercritical values 
all values satisfying the expression: 


Eq. yields larger values for than Eq. 32. The experimental values (1.25 
4.0) all are situated between our theoretical curves and 3.5, 
and thus, are below the curve (see Fig. 8). looks Eq. yields 
mean value all measurements both for stable and unstable flow conditions. 
Since almost all experimental values yield larger values for than 
those given Eq. one would question Tracy and Lester’s conclusion that 
contrary, that the measured values are reasonably good accordance with 
the computed values. 


Brock (2) did not find any influence the Froude number the friction 
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Brock’s experimental 


} 


—— Brock’s theoretical 
—— Authors’ theoretical 


(Betta) 


108 


FIG. Obtained Brock Comparison with Computed Values 


factor either. All measured values are situated the curve: 


which again yields values excess those for subcritical flow smooth 
channels (see Eq. 32). Brock mentioned (2) that error 0.5 (0.02 in.) 
All experimental values Brock are situated below the curve 
(see Fig. 9), and taking into account the error +4%, seen that they 
are good agreement with the values obtained from Eq. 42. The fact that 
Brock had only seven experimental results small interval Froude numbers 
(1.19 2.88) did not allow the construction curves for different values 
Moreover, seen from Fig. that the values for the two largest 
values (2.88 and 2.65) are distinctly larger than the values obtained from 
Eq. 51, which may not fortuitous. For the measurements 1.32 and 
1.50 2.2 10° and 8.4 10°, respectively) follows from Fig. that 
F.. These two values are the curve within 5%. 

the Hydraulics Laboratory the Ghent State University, Ghent, Belgium, 
some measurements (20) were carried out flume made acrylic sheet 
(40 ft) length and 0.24 (0.94 in.) width (1). The waterdepth was 
measured with point gage (three measurements each three cross sections). 
The discharge was measured means electromagnetic flowmeter. The 
relative error was: 
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0.01, one obtains Af/f from which results error 
2%-4% 

Measurements stable flow situations F.) agreed very well with Eq. 
(1). Some results unstable flow conditions are shown Figs. and 
11. These results are qualitatively good agreement with the theoretical results 
far the increase with concerned. Quantitatively, the measured 
values are approx 3%-4% smaller than the computed values, which 
difference the order the error the experimental result. 


7 

= 2.1 


—--- Van Heste’s experimental 


—— Authors’ theoretical 


5.103 alt 05 


FIG. 10.—Experimental Results Van Heste Constant Froude Number, 
Comparison with Computed Values 


Fig. that especially for large channels, the curves are very flat over large 
interval numbers. This might explain the fact that some experimenters 
working rather wide channel and restricted interval Reynolds numbers 
(2,19) could not recognize the effect the Froude number the friction 
factor, which was, the conditions their experiments, the order the 
error the measurements. 

Both Tracy and Lester (19) and Brock (2) worked with rather large channel 
(see Table and interval Reynolds numbers which almost 
independent (see Fig. 12). This might explain why they could not recognize 
the influence the Froude number the value but instead reported 
larger values for than given the Blasius Thijsse formula. For example, 
for 10° Brock found 0.0202, and Tracy and Lester found 0.0197, 
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sin@ = 0.015 
= 0.031 
“sine 0.039 


Vor Heste's experimenta 
—— Authors theoretical 


(B= 0.26m) 


—— 


FIG. Results Van Heste Constant Slope Angle, Compari- 
son with Computed Values 


FIG. 12.—Relative Increase, Function Reynolds Number, Channel 
Width, and Slope Angle 
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while Eq. yields 0.0182, which means increase 11% and 8%, 
respectively. These figures are good agreement with the results presented 
Fig. since reported channel slopes were less than 1°. Wakhlu, although 
working large channel well, did find influence because 
worked very large interval numbers (see Table 5). All other authors, 
the other hand, working narrower channels, found influence 


Using the White-Colebrook friction formula, the critical Froude number for 
turbulent flow becoming unstable can found function the shape 
the velocity profile, the Reynolds number and the relative channel roughness, 
the wave length the initial disturbance wave, and the channel width. For 
given Reynolds number (discharge), flow instability more likely occur 
wide rather than narrow channels. When scale models are built constant 
channels for unstable turbulent flow, might occur that the flow the 
prototype unstable while the flow the model remains stable. 

Introducing the critical Froude number the formula Rouse, the influence 
the Froude number the value the friction factor for unstable turbulent 
flow both smooth and rough channels can calculated. The value 
increases with increasing Froude number (channel slope). Moreover, found 
that the coefficients the formula relating and are functions the 
channel width. 


For large channels, the relative increase with respect results obtained 


with the White-Colebrook formula, becomes almost independent the Reynolds 
number, provided that This fact explains why some authors did 
not appreciate the influence although other authors, working narrow 
channels and for low numbers, did indeed report increase with 
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The following symbols are used this paper: 


cross-sectional area; 

channel width; 

dimensionless wave number 
Chézy roughness factor; 
dimensionless wave celerity; 

critical dimensionless wave celerity; 


448 
2 


UNSTABLE TURBULENT FLOW 


constants, the value which depends upon F(z); 
wave celerity; 
Froude number; 
critical Froude number; 
shape function the velocity profile; 
Darcy-Weissback friction factor; 
friction factor uniform flow, according White-Cole- 
brook; 
acceleration gravity; 
waterdepth; 
equivalent relative channel roughness; 
Manning roughness factor; 
wetted perimeter; 
pressure point (x,y), time 
discharge; 
hydraulic radius; 
Reynolds number; 
channel slope; 
time; 
mean velocity over cross section; 
velocity components the direction and direction; 
rectangular Cartesian coordinates; 
velocity distribution coefficient; 
wave number; 
roughness concentration; 
infinitesimal small amplitude the disturbance wave; 
slope angle the channel; 
dimensionless wave length 
wave length; 
kinematic viscosity the liquid; 
mass density the fluid; 
bottom shear stress; 
Reynolds stresses; and 
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SEPARATED FLOW THROUGH LARGE 
SUDDEN EXPANSIONS 


Prabhulal 


INTRODUCTION 


Most flows through sudden expansions closed conduits are turbulent with 
flow separation and recirculation. the expansions are axisymmetric (circular 
pipes), then recirculating regions and flow patterns through the expansions are 
symmetric. This has been observed and reported Kalinske (6), Chaturvedi 
(2), and Hue and Hebrard (5). When the expansions are rectangular ducts, 
depending upon whether the expansion ratio less than greater than 1.5, 
the flow patterns and recirculating regions have been found symmetric 
and steady, unsymmetric and unsteady, respectively. This has been reported 
Abbott and Kline (1), Filletti and Kays (4), and Durst, al. (3). Using 
analytical and experimental findings, the writer (7) has also reported that the 
flow patterns through rectangular ducts with expansion ratio 1.25 are 
symmetric and steady, whereas those with expansion ratio 2.5 are 
unsymmetric and unsteady. this paper, experimental findings for rectangular 
ducts with expansion ratios 2.0 and 3.0 are presented. These results will 
useful for predicting more accurately the behavior flows through large, 
sudden rectangular expansions. 


ConriGurRATION AND INSTRUMENTATION 


Problem Defined.—Fig. shows sudden rectangular conduit expansion 
ratio B/b, which and widths the rectangular conduit before and 
after the expansion, respectively; and The flow, after leaving 
the plane the expansion, develops two unequal cavities (stall regions) 
either side the axis symmetry, shown means full and dotted 
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lines Fig. The flow patterns which the reattachment point the lower 
boundary nearer (at the plane expansion shall referred 
flow type and that which the point reattachment the lower boundary 
away (at shall referred flow type II. Flow types and may 
get interchanged, giving rise unsteady behavior the flow. Between 
the primary flow the axial direction, and the secondary recirculating flow 
the stall regions, there exists free turbulent shear layer which highly 
sensitive disturbances, and causes these alternate types flow patterns. 
Experimental Setup.—A schematic representation the experimental arrange- 
ment shown Fig. 5-hp air blower was used supply air the 


Plane Sudden expansion 


- 
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FIG. 1.—Definition Sketch 
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FIG. Diagram for Experimental Setup; Expansion Ratio 2.0; All Dimen- 
sions inches (centimeters) 


inlet conduit with 2.5 cross section and 2.4-m length. Care 
was taken see that the air supplied was dust free and the test assembly 
free any mechanical vibrations. Flow the inlet conduit was fully developed, 
and the level initial turbulence was kept between 1%-2% the mean axial 
flow velocity, the inlet conduit. The inlet conduit was made join 
the expansion piece means detachable male-female joint, and the width 
the inlet and expansion conduit were the same (namely, 2.5 cm). This 
arrangement facilitated changes the expansion ratio any desired value. 
Four different expansion ratios (B/b 1.25, 2.0, 2.5, and 3.0) were used, 
and flow patterns were obtained for the Reynolds number range 0.5 10°-1.0 
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SHORT STALL 
REGION 


LEGEND 
MEAN VELOCITY 
MEAN PRESSURE 
PTS. OF MAX VELOCITY 
PTS. OF REATIACHMENT 


eeeeee MEAN VELOCITY FOR 
FLOW TYPE I 


FIG. 3.—Measured Values Mean Pressure and Mean Velocity; Expansion Ratio 
2.0; 10° (Flow Type and 
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FIG. 4.—Measured Values Mean Pressure and Mean Velocity; Expansion Ratio 
3.0; 10° (Flow Type 
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FIG. 5.—Measured Values Mean Pressure and Mean Velocity; Expansion Ratio 
3.0; 10° (Flow Type 
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10°. Experimental findings for expansion ratios 2.0 and 3.0 Reynolds 
number 10° are presented here. The Reynolds number was calculated 
which four times the hydraulic radius the conduit before the expansion. 
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FIG. 7.—Measured Values Longitudinal Turbulence Intensity; Expansion Ratio 3.0; 
10° 


Instrumentation.—A combined, static, and dynamic-type pitot tube was used 
measure the mean velocity, pressure, and wall shear the conduit downstream 
the expansion. Turbulence intensities the longitudinal direction, 
and the transverse direction, Vv’ and the Reynolds shear stress, 
were measured Dansk Industri Syndikat Associates (of Denmark) 
hot wire anemometer constant temperature type. The hot wire probe and 
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the bridge circuit measure cold resistance probe were developed the 
laboratory Bombay. The sensing element the probe was copper- 
coated tungsten filament diam, and was soft soldered the needle 
supports the probe. The central length the filament was etched 


RaTiO OF 


0-05 0-10 0-15,0 0-05 


Uo 


FIG. 


\ 


FIG. Values Transverse Turbulence Intensity; Expansion Ratio 2.0; 


out, making use nitric acid, give cold resistance the sensing element 
around ohms. 
straight wire probe was used measure and single, 45° 
inclined wire, which could rotated through 180°, was used measure 
The Reynolds shear stress, was calculated from the 
instantaneous readings and v’. 
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FIG. Values Transverse Turbulence Intensity; Expansion Ratio 3.0; 


FIG. 


456 HY4 


HY4 LARGE SUDDEN EXPANSIONS 
FINDINGS 


Mean Velocities and Pressures.—Figs. and show measurements the 
mean velocity profile and static pressure distributions various locations 
downstream the plane expansion for expansion ratios 2.0 and 3.0, 
respectively. The degree asymmetry larger for the expansion ratio 3.0 
than for the expansion ratio 2.0. Maximum mean velocities would expected 
occur the center line the conduit all locations, but the measurements 
show that the line connecting points maximum mean velocity deviate from 
the conduit center line. The maximum deviation from the center line for expansion 
ratios 2.0 and 3.0 0.2 and 0.466 respectively. Flow patterns also 
show unequal reattachment lengths. For expansion ratio 2.0, the reattach- 
ment lengths are found equal and whereas for the expansion 
ratio 3.0, they equal 4.5 and 15.0 The flow becomes symmetric and 
stable (maximum value mean velocity the duct center line) x/T 


FIG. 10.—Measured Values Turbulent Shear; Expansion Ratio 2.0; 10° 


20.0 for expansion ratio 2.0, and for expansion ratio 
3.0. Figs. and show that flow types and are mirror images. 

Mean static pressures distributions also show greater asymmetry for larger 
expansion ratios (similar the mean velocity patterns), and also become 
symmetric (downstream the plane expansion) the same locations 
the mean velocity distributions. 

Turbulent Characteristics.—Figs. and display the results the turbulence 
intensity the longitudinal direction, and Figs. and show the results 
the turbulence intensity the transverse direction. Measured values the 
turbulent shear stress are given Figs. and 11. 

Turbulence intensities and shear stress show higher degree nonuniformity 
for the expansion ratio 3.0 than for the expansion ratio 2.0. Peak values 
deveiop along the free turbulent shear layer. Secondary peaks also occur 
the mideddy pocket region nearer the boundary. The turbulence grows more 
rapidly and decays more quickly for expansion ratio 3.0 than for 
expansion ratio 2.0. Fig. shows that flow types and are mirror images. 


= é 
28 
RATIO OF x/T 12 } 0-7 
i 10 | $. 
| % \ | | 6 
"0 
150 100 50 0,200 150 100 50 0 
++ x10 10° 
Ug’ Uo? ve 
7 


APRIL 1981 


Uo* 


| RATIO OF x/T 
| 


120 


uv 


FIG. 11.—Continued 


458 HY4 
| | 6 
160 120 60 20 0 200 160 60 40 
FIG. 11.—Measured Values Turbulent Expansion Ratio 3.0; 10° 
91-0 | 
T=10 
| 
1-0 


HY4 LARGE SUDDEN EXPANSIONS 


The experimental data given this paper provides useful information regarding 
the behavior the mean flow pattern and the turbulent characteristics for 
two-dimensional flows through large, sudden expansions. These results comple- 
ment and extend the previous results (7), and show that the expansion ratio 
affects all the important parameters the flow field. More specifically, with 
increase the expansion ratio, the reattachment lengths are found decrease 
for the short stall region and increase for the long stall region. For larger 
expansion ratios, the flow patterns and flow characteristics the two types 
flow (types and II) are found mirror images each other. The 
rate development the peak values turbulence and the rate its decay 
are found higher for larger expansion ratios. The Reynolds number does 
not influence the flow patterns for the range studied (0.5 10°-1 10°). 
The flow pattern becomes symmetric and stable sooner, and the rate mixing 
higher for larger expansion ratios. 


TABLE Expansion Ratio Turbulent Characteristics 


Percentage Values Level Turbulence 
Longitudinal Direction 
Expansion ratio and Maximum Uniform 
Reynolds number value value 


(1) 


With increase the expansion ratio, mean flow patterns become more 
unsymmetric and unsteady (see Table Ref. 7), whereas the degree level 
turbulence remains almost unaltered except that the peak values develop 
earlier and decay faster than those with low expansion ratios (see Table 1). 

The results experimental findings for large expansion ratios will useful 
designing those engineering devices which high rates natural mixing 
are desired, e.g., chlorine and alum injecting plants water works, color mixing 
plants chemical processes, etc. Present experimental findings will also 
useful constructing suitable mathematical model for predicting the behavior 
large sudden expansions. 
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The following symbols are used this paper: 


B,b conduit width before and after expansion, respectively; 
Static pressure; 
Reynolds number (given which four 
times the hydraulic mean depth the inlet); 
step height given b); 
mean axial velocity the inlet; 
turbulent shear stress; 


root mean square values and respectively; 
coefficient dynamic viscosity; 

mass density; 

denotes temporal mean; and 

denotes instantaneous deviation from the temporal mean. 
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SIMULATING TRANSPORT 
UNCONFINED AQUIFERS 


Varuttamadhira and Raymond Volker,’ ASCE 


more water extracted from surface and underground sources satisfy 
increasing water demands, more polluted water returned which may contaminate 
those sources. Shallow unconfined aquifers are particularly susceptible 
pollution from conservative contaminants when little natural treatment 
afforded overlying strata. Problems both steady and transient flow 
unconfined aquifers are often solved considering only the saturated domain 
below the water table (4,17). Alternatively, the unsaturated zone may also 
included (13,16). Mass transport unconfined aquifers has also been treated 
both methods (14,18,19,20,21,22). The approach, using only the saturated 
region, usually less expensive because the problems nonlinearity the 
governing equations are less formidable even though the solution domain may 
nonstationary. 

This paper presents solution the mass transport problem for unconfined 
aquifers. The numerical scheme employed the conventional weighted residual 
finite element method (24) with special provisions allow for the moving water 
table. The finite element mesh changed iteratively locate the new water 
table each time step when solving the flow equation. Thus, the solution 
the mass transport equation, the nodes the mesh may have new positions 
successive time intervals. The solution technique presented eliminates the 
need for interpolation solute concentration values accomodate this movement 
the mesh. 
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The continuity condition ground-water flow can tensorially expressed 


which velocity the direction x,. (Note the use the Einstein 
summation convention.) Using Darcy’s law, Eq. can written: 


which hydraulic conductivity tensor; and piezometric head. 

The problem studied that surface basin recharging unconfined aquifer, 
depicted Fig. The basin width 2L, and the aquifer width between 
prescribed hydrostatic pressure boundaries 2B. The basin and aquifer are 
assumed sufficiently long for the flow two-dimensional. With the 
aquifer water level initially above the base, recharge rate causes 


FIG. 1.—Definition Sketch 


the water table rise with time new height which dependent the 
distance from the center line the basin. Since the field symmetric, only 
half included the solution and the origin coordinates taken 
the base the aquifer below the center line the recharge basin. The 


flow governed Eq. which subjected the following boundary and 
initial conditions: 


h=a, for t=0, x,€(0,B), 


for x,=0 


onthe free surface, where the height the free 


HY4 

(3): 
f 


HY4 SIMULATING MASS TRANSPORT 
h=x, for x,=B, 


with R=P, for and for 


which initial saturated depth; aquifer half width; basin half 
0); direction cosine outward unit vector normal the free surface; 
time; infiltration rate; and infiltration rate from the basin. 


The hydrodynamic dispersion equation for incompressible flow given 
Bear (3) as: 


(see Ref. 7); and hydrodynamic dispersion coefficient tensor. 
Eq. subjected the following conditions: 
—=0 for x,=0, and 


for 


for 


which initial moisture content the saturated zone; and 
concentration pollutant the unsaturated zone. Eq. 12, which has been 
found work well previous studies (21,22), allows convection solute 
where there velocity normal the boundary. The derivations Eqs. 
and are presented Appendix should stressed that the formulation 
presented Appendix based the assumption sharp free surface. 


Furthermore, Eq. can only macroscopically describe the mixing the free 


(9) 
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surface, because the assumption that the mass influx instantaneously mixed 
with the mass already the saturated zone. Although this assumption contradicts 
the experimental findings Warrick, al. (23), which indicate that the advance 
the saturated front into the unsaturated zone precedes that the solute 


front, the results obtained with the present analysis give reasonable agreement 
with the experiment. 


Finrre 
Flow Equation.—Assuming that the solution can represented by: 
NN 


which approximate solution; value state variable node 
element. 


Applying the weighted residual technique Eq. one obtains 


which area over which integration takes place; and weighting 


function. 
Using Green’s Theorem and 17, Eq. becomes: 


aN, 
A s 


Ox, 


which boundary over which integration takes place. 
Eq. substituted into Eq. give 


dh, 


which should noted that the free surface and therefore, 

Using the Galerkin technique, which N,, one obtains set 
simultaneous equations written matrix form as: 


e 


s 
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which represents the sum elemental contributions. 


Approximating the temporal derivative simple finite-difference approach, 


+ ) {h}'**! + (1 —r){F}‘ 


which time-stepping factor ranging from The writers (11) recommend 
value for good accuracy and stability. 

The system simultaneous equations solved iteratively until convergence 
reached, i.e., for all nodes the free surface, the maximum difference 
between the piezometric head and the coordinate the node must less 
than some stipulated tolerance. 

Hydrodynamic Dispersion Equation.—In order accommodate the expansion 
contraction the saturated domain, the solution sought through the finite 
element method with moving mesh. using such mesh necessary 
modify the transport equation take account the mesh movement. 

Let the equations transformation from curvilinear coordinate system 


rectangular cartesian coordinate system (where the coordinates point, 
are 


which curvilinear coordinates the system moving relative the 
rectangular system. Eq. can inverted give 


t), 


which valid provided that the Jacobian, does not vanish (2) where 
defined 


Mollowney (15) has shown that the transport equation for compressible 
fluid the frame reference 


— + = y 


can written the frame reference 


é 


which terms with asterisks are equivalent their counterparts the frame. 
Relationships between variables these two reference frames are 
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Therefore, from Eq. 34, the velocity the moving frame relative 
the stationary frame observed the frame. 

The Jacobian neither spatially nor temporally invariant and this gives rise 
difficulties accurately evaluating the frame reference. overcome 
this problem, assumed that each element geometrically time invariant 
each time step with geometry which the average those time 
and At. This equivalent assumption unit Jacobian for each 
element which implies that the motion isochoric. 

The preceding true the following transformation employed between 
tand 


which time elapsed from and assumed constant over 
the subdomain between and Ar. 
From one can verify that each subdomain 


conversely 


which Kronecker delta. 
From the assumption that constant each subdomain, and from the 
condition incompressible flow can established using Eq. that 


Using Eq. and the condition isochoric motion 1), Eq. can 
further reduced 


ac* ac* ac* 


For the present problem, the solution sought the plane, and 
approximated the form: 


NN 


Applying the weighted residual method Eq. one obtains 
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ac* ac* ac* 
Using Green’s theorem and Eq. 39, Eq. becomes 
aN, 
dt A 0&, 
P ae, i 


which yields system simultaneous equations written matrix form as: 


[MM] | [STT] {FF}‘**/> =0 


e A 
e s 


Using 30, 31, 32, 33, 34, and 36, can stated that each element 


From the preceding relationships one can see that coincides with 
and coincides with Their values calculated from the frame can therefore 
conveniently utilized the foregoing finite element formulation. 

The velocity, w,, each element approximated which Ax, 
the movement the centroid the element between time and 
obtained the solution the flow equation and resulting from the free-surface 
movement. 


Replacing the and terms 


Eq. one obtains 
[MM 1 t+(At/2) 


l MM t+(4t/2) 
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This equation uses time-stepping factor 0.5 which, for the case Galerkin’s 
method (N, W,), the stationary domain, has been shown Gray and 
Pinder (9) minimize smearing the front the expense numerical dispersion. 


Experiments for situation similar that shown Fig. were conducted 
sand box model verify the mathematical models. Because symmetry 
only one half considered here. 


Data Acquisition System 


Sand Box Model 
FIG. 2.—Layout Experimental Facilities 


Node number 


FIG. 3.—Mesh Layout (Nonstationary Element) 


matrix conductivity probes was embedded the sand box various 
locations monitor the movement the contaminated front. These probes 
were connected conductometer which was turn linked with minicomputer. 
Piezometric heads various locations along the sand box were recorded manually 
piezometric tubes and automatically, using pressure transducers which were 
also connected the minicomputer. This data acquisition system diagrammati- 
cally shown Fig. 

each experiment, solution was employed represent conservative 
pollutant and was introduced the upstream end the sand box 


POP-8/E PDP-10 
Mini Computer Main 
Computer 
Pressure 
Conductometer 
Water Recharge box 
Supply Information 
System 
— 
10 
0 q 0 
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L). The concentration the solution was the order 0.1% weight, 
that the change density and viscosity due the presence the salt 


solution was negligible. For more details, readers are referred the first writer 
(10). 


The accuracy the nonstationary element scheme was first tested using 
one-dimensional model hypothetical situation which the domain expands 
with time. shown Fig. there constant concentration source 
The right-hand edge the flow domain arbitrarily assumed move with 
velocity, 4u, and this turn causes the discretization inside the domain 
expand accommodate the new domain. The left-hand side the domain 
remains stationary all times, and the boundary condition that point 
C(o,t) C,. analytical solution this particular problem 


which erfc (x) complimentary error function argument, 

comparison between the analytical and numerical solutions made Fig. 
from which evident that the numerical solution agrees well with the 
analytical solution. 


Analytical 
° Finite Element 


12 24 36 


FIG. 4.—Concentration Profiles 


typical discretization the nonstationary two-dimensional domain used 
287 nodes rows columns) and 480 linear triangular elements with the 
horizontal node spacing increasing with distance from the recharge basin. Since 
the horizontal length aquifer fixed, the water table movement can 
simulated using horizontally-constrained elements with the nodes moving 
only vertically. 

Data from two experiments were selected verify the two-dimensional 
numerical models. shown Table the original saturated depth experiment 
larger than that experiment while the reverse true for infiltration 
rates. consequence, the rise the free surface relation the original 
saturated depth larger experiment 


0.0 
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The coefficients advective dispersion for the homogeneous sand medium 
were calculated first estimating longitudinal and lateral coefficients advective 
dispersion, and D,, respectively along the principal axes using the following 
expressions which were obtained from laboratory experiments (10): 


which the Reynolds number, given 


which local pore velocity along the instantaneous streamline; 
average grain-size diameter; and kinematic viscosity. The numerical coeffi- 


TABLE 1.—Parameter Values for Experiments and 


milli- milli- 
milli- per per 
meters meters meters second second 


(3) (4) (5) (6) (7) 


FIG. 5.—Comparison Measured and Calculated Free-Surface Levels for Experi- 
ment 


cients Eqs. 48a and are good agreement with those obtained other 
workers for similar laboratory experiments. The components, D,, were then 
found tensorially rotating and back the reference axes (12). 

Verification the flow model shown Fig. which free-surface 
(water table) levels various locations are compared. Results shown Fig. 
indicate that the flow model can accurately simulate the state variable, 
piezometric head, inside the flow domain. For further details the flow model, 
readers are referred Ref. 11. 


v 
(1) (2) (8) (9) 
0.1645 6.0 0.30 1.8 1.009 0.325 0.491 
0.308 6.0 0.30 1.8 0.445 0.303 0.491 
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Verification the transport model was made comparing the experimental 


and calculated breakthrough curves shown Figs. and Fig. results 
are drawn from two cases incorporating two different boundary conditions 


the free surface outside the recharge basin; namely: 


case2: 
equated zero since the moisture the unsaturated zone 


case 
that freshwater. Thus, application the boundary condition case 
x,/L = 14.99 


Legend 


Experiment 
Finite Element 

with no Dilution (Case u 
Finite Element 


~~ with Dilution (Case 2) / 


18 
t Po/ (8, ae) 


FIG. 6.—Comparison Experimental and Calculated Breakthrough Curves for Experi- 


ment x,/a, 1.70 


Experiment 


Finite Element 


with no Dilution } 


FIG. 7.—Comparison Experimental and Calculated Breakthrough Curves for Experi- 
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c/c, contours 


No Dilution at Free Surface 


_ With Dilution at Free Surface 


Experiment 


Finite Element 
with no Dilution 


FIG. Experimental and Calculated Breakthrough Curves for Experi- 
ment 


C/Ce contours 


No Dilution at Free Surface 


With Dilution at Free Surface 


FIG. Concentration Contours for Experiment 
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implies that dilution the pollutant the saturated zone takes place the 
vicinity the free surface. Concentration contours from the preceding two 
cases are exhibited Fig. which one can see that the inclusion the 
moisture the unsaturated zone results the retardation the front the 
vicinity the free surface. 

Figs. and show results from experiment inspection Fig. 
which shows breakthrough curves locations close the free surface reveals 
that the use the boundary condition Eq. causes the breakthrough curves 
lag behind those employing the boundary condition Eq. 50. Fig. 
x,/L 14.99, dilution the front due the inclusion moisture from 
the unsaturated zone can observed the experimental result, although not 
the extent predicted case Fig. the predicted piezometric heads 
are somewhat lower than those experimentally measured. This implies that outside 
the recharge basin (x, L), the calculated free surface also lower than 
its experimental counterpart, because the pressure distribution almost hy- 
drostatic that region (10). Calculated breakthrough curves Fig. resulting 
from Eq. 51, therefore appear suffer more dilution than the experimental 
curves. Also Fig. the calculated breakthrough curves not correspond 
well the experimental ones 0.39, which under the infiltration 
basin. The calculated breakthrough curve abrupt because the abrupt rise 
the calculated free surface passing the probe location; however during the 
experiment, the probe did register the concentration the NaCl infiltrate while 
above the free surface. 

The inclusion the moisture influx across the free surface does not affect 
breakthrough curves locations well inside the flow domain, and only results 
derived using Eq. the boundary condition are compared against the 
experimental breakthrough curves Fig. 

Results from the finite element transport models are compared with those 
from experiment Fig. this case, breakthrough curves derived using 
either Eq. Eq. boundary condition are almost identical. Only those 
using Eq. are displayed Fig. From the concentration contours shown 
Fig. 10, apparent that the moisture influx the free surface does 
not appreciably affect the mass distribution inside the saturated domain. 


ANALysis 


Crank and Gupta (6) proposed moving grid system solve the one-dimensional 
diffusion problem involving moving boundary. Their method employs 
finite-difference grid which moves with the boundary and utilizes cubic spline 
Lagrangian interpolation schemes determine values concentration from 
the previous timestep the new locations. Al-Niami and Rushton (1) also proposed 
similar moving grid technique using linear interpolation. Such forms 
interpolation are not practical multidimensional problems, especially those 
using finite element mesh with geometrically irregular elements. The numerical 
scheme proposed this study surpasses the moving grid system this respect, 
because does not require any kind interpolation between timesteps. 

Figs, and indicates that there exist some discrepancies 
between the models and experimental results. There are several possible explana- 
tions for this, including experimental errors such pockets nonhomogeneity 
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the sand, possible fingering instability the salt front, and inherent numerical 
errors the numerical models. The differences between the numerical and 
experimental results are relatively small, and could accounted for the 
aforementioned factors. may therefore inferred that the flow and transport 
models are accurate simulating the movement mass the saturated zone. 

From Figs. and 10, evident that the magnitude the area which 
the mass distribution affected the moisture flux across the free-surface 
boundary dependent upon the rate rise the free surface which, 
turn, governed the infiltration rate. Thus, for large rise the free 
surface, lack knowledge the mass the unsaturated zone may result 
appreciable error predicting the mass movement the neighborhood 
the free surface. 


this study, the weighted residual technique, conjunction with the finite 
element method, was utilized formulate numerical models flow and transport 
nonstationary saturated ground-water zone. Through the formulation present- 
ed, elements are allowed expand and contract required accommodate 
the nonstationary saturated domain without the necessity for interpolation 
solute concentration values when proceeding from one timestep the next. 

Comparisons between the numerical results and those from limited number 
experiments indicate that the models can accurately simulate the flow 
water and the movement pollutants the saturated zone homogeneous 
isotropic aquifer. Even though the models are primarily intended for the saturated 
zone, they can also approximate the mass transfer result the movement 
the saturated-unsaturated interface (free surface). The models can thus 
linked with independent unsaturated flow and transport models areas where 
the flow the unsaturated zone significant. This should provide more 
economical alternative the conventional saturated-unsaturated scheme. 


Considering elemental volume the free surface depicted Fig. 
during time At, the principle conservation mass written: mass 
retained net contribution mass from which 


which concentration contaminant the free surface. the experimental 
work, the concentration actually measured the recharge basin and not 
the water table. This should not lead significant errors because the dispersion 
during movement down the water table not large, and once the flow 
fully established that zone, the solute concentration the water table will 
effectively that the recharge source. 
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The preceding equation implies that the invading fluid instantaneously mixed 
with the resident fluid the free surface. The basic concept pore space 
division shown Fig. from which follows that 


This assumes that the stagnant dead-end portion where fluid usually 
considered immobile (5,8) negligible, that the Darcy porosity (or 
interconnected porosity) approximately equal the total porosity. 

The first term the right-hand side Eq. denotes the influx due 
infiltration, the second term the flux from the saturated domain, and the last 
the contribution mass already occupying the pore space. 


(b) moisture 
content 


-p, Kij ah 


FIG. 11.—Free-Surface Definition: (a) Elementary Volume; (b) Hypothetical Moisture 
Profile Free Surface 


assumed that the moisture the portion has been replaced with 
tracer concentration 


Dividing Eq. through and using Eq. Eq. becomes 


which 14. 


Similarly, where and tracer concentration occupies the water 
the unsaturated zone which has moisture content one can easily 
verify that Eq. results. 

For the case receding free surface outside the basin, Eq. 
employed because the concentration immediately inside the free surface 
equal that immediately outside. 
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The following symbols are used this paper: 


n 


initial saturated depth; 

half width aquifer; 

concentration; 

concentration the unsaturated zone; 
concentration the inflow boundary; 
hydrodynamic dispersion coefficient tensor; 
longitudinal convective dispersion coefficient; 
lateral convective dispersion coefficient; 

mean grain-size diameter; 

free-surface elevation; 

piezometric head; 

Jacobian; 

hydraulic conductivity tensor; 

half width basin; 

component unit outward normal the boundary; 
basis function node 

infiltration rate; 

ith component superficial velocity 
recharge rate; 

time; 

absolute pore velocity; 

ith component pore velocity; 

weighting function node 

velocity the moving frame relative the stationary frame observed 
the stationary frame; 

cartesian coordinate the stationary frame; 

time increment; 

Kronecker delta; 

effective porosity; 

initial moisture content; 

Darcy porosity; 

kinematic viscosity; 

moving coordinates; and 

W,. 
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INTRODUCTION 


Modeling unsteady shallow water flow porous media has many engineering 
applications, including the investigation overland flow, infiltration, surface 
irrigation, and seepage from channels. Theoretically, the process described 
conjunctive surface-subsurface flow system. The subsurface flow region 
may contain both saturated and unsaturated zones. The surface and subsurface 
flow components are dynamic equilibrium and they interact through the process 
infiltration. 

Surface flow was studied separate component the past. number 
researchers including Wooding (30), Chen and Chow (4), and Chow and 
Ben-Zvi (5) proposed various methods simulate surface runoff which 
infiltration was either disregarded treated somewhat casually. Also, 
mathematical models developed for surface irrigation including those Chen 
(3), Hart al. (10), and Sakkas and Strelkoff (20), infiltration rates were 
approximately estimated. Meanwhile, most subsurface flow studies were restrict- 
either unsaturated zone saturated ground-water flow. Whisler and 
Klute (29), Staple (24), Remson al. (18), and Rubin (19) studied vertical 
moisture flow unsaturated soil column. 

Verma and Brutsaert (28), Prickett and Lonnquist (16), and Singh (22) examined 
saturated ground-water flow disregarding the zone aeration. Saturated and 
unsaturated zones were combined Freeze (7), Hornberger and Remson (11), 
Abbott al. (1), and Vauclin al. (27). Reeder al. (17) examined the effect 
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variation water depth infiltration treating the surface water depth 
prescribed function time. conjunctive surface-subsurface flow model 
was first developed Smith and Woolhiser (23) remarkable study 
overland flow pervious surfaces. However, they simplified the problem 
using kinematic wave approximation represent the surface flow and 
considering the subsurface flow with only vertical movement water through 
unsaturated soil column bounded fixed water table the bottom. 

more sophisticated conjunctive surface-subsurface flow model developed 
the present study. Surface flow described set one-dimensional 
dynamic wave equations. Subsurface flow assumed two-dimensional 
with potential gradients the vertical well the surface flow direction. 
Both saturated and unsaturated zones are considered. All components the 
system are dynamic equilibrium, satisfying the principles continuity and 
momentum. brief description this model given this paper. Further 
details can found Ref. 


Surrace 


the present study, surface runoff represented the Saint-Venant equations 
simplified for wide channels 


+ 


flow depth normal time; time rate lateral inflow; time 
rate infiltration; gravitational acceleration; bottom slope; and 
friction slope. The assumptions and the rigorous derivations the 
equations were reported Strelkoff (25) and Yen (32) and will not repeated 
here. 


The friction slope computed using the Darcy-Weisbach formula 


The frictional resistance coefficient evaluated from approximate Moody 
diagram utilized watershed model proposed Kareliotis and Chow (13). 
Depending upon the instantaneous state the flow, evaluated using Eqs. 
4-6 respectively, for laminar, transitional, and fully turbulent flow 


y? 
0.223 
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1.74 


Eqs. 4-6, the Reynolds number, equal the product divided 
the kinematic viscosity; length measure surface roughness; and 
the equation 


which rainfall intensity in./h. The effect the pervious ground 
surface leaking boundary negligible for the range infiltration 
velocity surface velocity ratio investigated the present study (31). 

four-point implicit finite difference scheme employed for the solution 
the surface runoff equations. shown Fig. the flow domain divided 
into number flow reaches. Similarly, the time variable also discretized 


FIG. 1.—Schematic Diagram Surface Flow and Computational Grid 


and solutions discrete time intervals are sought. Accordingly, the Saint-Venant 
equations (Eqs. and are written finite difference quotients 
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(yrs)? 


+ + re (l +y?) —y;) 

+ 

which fixed computation time interval; the projection horizontal 
plane the length the flow reach which centered the ith grid station; 
the average time rates rainfall and infiltration, respectively, 
along the flow reach between ith and 1)th grid stations; and weighting 
factor. The superscript refers time levels computation. The subscript 
refers flow sections along the flow length. 

using Eqs. The lateral inflow known from rainfall surface 
water supply data. The rate infiltration, determined from the coupling 
surface and subsurface flow equations will considered later. All the 
other terms Eqs. and are known either from initial conditions from 
the previous time-step computations. 

Initially, the flow depth and velocity must known all the flow sections 
the initial condition. the case overland flow developing initially 
dry surface, very thin water film assumed ponded before the flow 
initiated. This assumption necessary overcome the numerical singularity 
the solution procedure. The results are not affected the assumed thickness 
this very thin film (2). One boundary condition each the upstream 
and downstream ends the surface domain needed for subcritical flow. 
The boundary conditions may prescribed discharge stage hydrograph 
discharge-stage relationship. 

applying Eqs. and each the rectangular grids along 
the surface flow length and combining these with two equations representing 
the boundary conditions, set simultaneous nonlinear algebraic equations 
(2N) unknowns obtained. The unknown quantities are and for 
generalized Newton iteration technique. The coefficient matrix obtained 
the iteration process has banded structure, and this property used great 
advantage devising fast solution technique (2). 


Subsurface flow modeled two-dimensional motion single-phase 
incompressible fluid nondeformable porous medium and represented 
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the equation 


relative permeability; piezometric head; capillary pressure head; 
soil porosity; degree saturation; and time. The derivation 
and physical interpretation Eq. was given elsewhere, e.g., Ref. and 
not repeated here. 

applied saturated groundwater flow well unsaturated soil 
water flow. the zone saturation the positive pressure head and 
unity. neglecting the entrapped gases the soil, for saturated flow 
also set equal unity and aS/aP 

Hydraulic relationships among and are needed order employ 
Eq. the zone aeration. Unfortunately, unique P-S-K, relationships 
cannot obtained since soils exhibit different behavior during wetting and 
drying processes. The present model partly accounts for this property soils 
known hysteresis allowing two different P-S-K, relationships: one for 
each the wetting and drying processes. 

finite difference scheme based the successive line over-relaxation 
method (SLOR) employed for the solution subsurface flow equation. 
shown Fig. the entire porous medium being analyzed divided into 
rectangular volume elements. Assuming soil properties remain constant within 
each element, Eq. written finite difference form for interior element 


Ax, 


( i,j +0.5 as AZ, +Az,, 
Az -1 + Az, 


n+ n+ n+ Az,, 


which volume the rectangular element enclosing the node (i, which 
the centroid the element; and the side lengths the element; 
time stage computation; and computational time increment. 

the number elements vertical strip shown Fig. The coefficients 
K,, and are prescribed functions All other terms 


ij 
are known either from specified initial condition from previous time-step 


= = 
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computations. Any antecedent moisture distribution can adopted the initial 
condition the subsurface flow. 

all external boundaries, Eq. modified incorporate the boundary 
conditions. Specified-flux specified-head conditions, which may time 
varying, constitute the boundary conditions the bottom and end boundaries 
the flow domain. The top boundary condition determined from the surface 
flow depth. the case initially dry surface, specified mass-flux boundary 
employed prior surface saturation. After thin layer soil becomes 
saturated the land surface, ponding overland flow initiated, and the 
top boundary condition changed from the specified mass-flux specified 
head condition. The pressure head each boundary element determined 
the surface water depth above it. 


FIG. Mesh for Subsurface Flow Solution 


nodes together with (2M boundary equations form set 
simultaneous nonlinear algebraic equations with unknowns, 
nonlinear equations obtained using successive line over-relaxation 
(SLOR) technique oriented the direction. Details the solution procedure 
was given Akan (2). 


Surface and subsurface flow components are interrelated common pressure 
head and exchange water the ground surface shown schematically 
Fig. The inter-relation between the two components depends upon the 
instantaneous flow conditions, and for transient flow changes with time. 
the mathematical sense, the top boundary condition for the subsurface flow 
domain determined the surface flow conditions. Conversely, infiltration, 
lateral outflow term Eq. controlled the soil water conditions. 
Therefore, solution can obtained only the finite difference equations 
representing the two components are coupled each time step. quasi-implicit 
coupling procedure adopted and described follows: 
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Using the subsurface flow solutions obtained specified the previous 
time stage, the rate infiltration each one the subsurface volume 
element the ground surface computed from 


aP él 


i,1+0.5 4,2 Az,,+A4z,, 
n+ n+ n+ Az —i1 + Az,, 

5x r7i+0.5,1 it +Ax,, 

Knowing the rate infiltration each surface flow cross section, the 
surface flow equations are solved. 

Using the surface flow solutions, the piezometric head each computational 
node lying the upper boundary the subsurface domain determined from 
the surface flow depth shown Fig. Employing these values the 
top boundary condition the current time step, the subsurface flow computed. 


Woter Surfoce 


Ground Surface 


infiltration 


° 
Pressure Head 


FIG. Diagram Pressure Head for Surface and Subsurface Flow 


With the moisture distribution computed step the rate infiltration 
re-estimated from Eq. 12, and the surface flow recomputed using the 
new estimates for infiltration. 


Steps 2-4 are repeated until the surface and subsurface flow solutions change 
within specified acceptable error for two consecutive iterations. However, 
usually one cycle iteration sufficient since the response the subsurface 
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flow the variation surface flow depth much slower than the response 
the surface flow changes the rate infiltration. 


The important assumptions involved the development the proposed model 
are listed below. 

For the surface flow, the flow gradually varied with respect space, 
accordingly: 


The piezometric pressure distributed uniformly over flow cross section. 

The velocity distributed uniformly over flow cross section. 

The effects boundary friction and turbulence surface flow can 
described 3-6. 


For the subsurface flow: 


Darcy’s Law applicable transient subsurface flow. 

The flow air phase the zone aeration negligible, and there 
effect air compression the liquid flow. 

The porous medium nondeformable, and the fluid incompressible and 
homogeneous. 

The osmotic pressure, adsorption, and thermal potential gradients are 
negligible. Alternatively, adsorption and osmotic pressure can treated 
part the capillary potential. 

The P-K,-S relationships not change with time. 


Naturally, these assumptions would restrict the usefulness the conjunctive 
model developed. However, from practical viewpoint, despite the assumptions 
the model can still used solve wide variety surface-subsurface flow 
problems. The unsteady, gradually varied surface flow can laminar, transitional, 
fully turbulent, uniform nonuniform, and subcritical supercritical provided 
roll waves similar complications developed. The transient subsurface flow 
can one- two-dimensional. The porous medium can homogeneous 
nonhomogeneous, isotropic anisotropic, and saturated unsaturated, both. 


order verify the proposed model, first the surface and subsurface flow 
components are tested separately with available experimental data and analytical 
solutions. Then, using set published experimental data, the model 
verified applied conjunctive surface-subsurface flow. 

Izzard’s (12) experimental results are adopted verify the surface flow 
component. Run 138 his experiment, the runoff from simulated uniform 
rainfall over the entire area plot was recorded. The plot was (1.8 
wide and (21.9 long and had bottom slope 0.001. The surface 
the plot was covered with crushed slate roofing paper which represented 
time-variable intensity the simulated rainfall, the observed runoff, and the 
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computed runoff are shown Fig. The sudden jump the runoff hydrograph 
when the rainfall terminates apparently due sudden decrease the 
flow resistance. can inspected visually Fig. the agreement between 
numerical and experimental results good. 
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0 
Time (min) 


FIG. 4.—Verification Surface Flow Prediction with Experimental Data 
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FIG. 5.—Verification Subsurface Flow with Analytical Solution 


order verify the subsurface flow component the proposed model, 
Philip’s (15) analytical solution for vertical flow water initially unsaturated 
soil, Luthin and Holmes’ (14) and Hantush’s (9) analytical solutions for horizontal 
saturated flow, and set experimental data reported Vachaud and Thony 
(26) are adopted. 


the idealized case vertical percolation water through deep homoge- 
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neous soil column, the initial saturation the soil assumed uniform and 
equal 0.20 throughout the entire depth. The surface soil becomes fully saturated 
zero time and remains throughout the entire process. Sufficient amount 
water supplied continuously the soil surface under zero pressure head. 
The idealized hydraulic properties the soil are K,, 6.25 mm/s, 
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FIG. 6.—Verification Subsurface Flow with Luthin and Holmes’ and 
Analytical Solutions 
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profiles obtained for this idealized situation using the proposed model and Philip’s 
analytical method are good agreement. 

another hypothetical situation, the subsurface flow component the 
proposed model applied the horizontal movement saturated ground-water 
flow. The initial elevation the water table above impervious base 
The variation the water table elevation examined the piezometric 
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FIG. 8.—Comparison Subsurface Flow Solution with Experimental Data and 
Another Numerical Solution 


Holmes’ (14) and Hantush’s (9) analytical solutions, the aquifer characterized 
storage coefficient 0.45 and permeability 0.36 m/h. The same 
aquifer reproduced the numerical scheme with 0.36 m/h, 
and the product aS/aP 0.05 can inspected Fig. the 
numerical solution good agreement with the analytical solutions. 

Vachaud and Thony’s (26) experimental study one-dimensional vertical 
movement water percolating through homogeneous unsaturated soil column, 
the soil moisture was initially static equilibrium. Throughout the process, 


° 
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the pressure head the soil surface was maintained —0.12 outflow 
occurred from the bottom the soil except the later part the experiment 
which not considered here. The hydraulic properties the soil measured 
and used the numerical simulation are given Fig. order avoid 
possible time lag between the initial times the experimental and the 
mathematically simulated conditions, the numerical solutions are started 
100 sec with the corresponding moisture distribution determined experimen- 
tally. comparison the saturation profiles obtained from the proposed model 
and from another numerical solution method Giesel al. (8) with the 
experimental results shown Fig. The agreement satisfactory. 
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FIG. Computed and Observed Saturation Profiles for Conjunctive 
Model 


conjunctive surface-subsurface flow system, the presently developed 
model tested with the experimental data Smith and Woolhiser (23). 
laboratory watershed consisting soil flume and rainfall simulator was 
employed the experiment. The soil flume was in. (51 mm) width, 
(1.22 depth, and (12.2 length. The bottom the flume 
was impervious, but seepage was allowed both ends. The thickness the 
three zones soil placed layers the flume are shown Fig. The 
hydraulic properties the soil layers are given Fig. 10. The fluid used 
was light oil with kinematic viscosity 0.003 in./sec (1.94 
The intensity the nozzle generated rainfall was 9.9 (250 mm/h) 


Soil Zone 


HY4 SHALLOW WATER FLOW 491 


over min. surface roughness 0.05 in. (1.3 mm) used 
the numerical simulation. The numerical and experimental results the soil 
moisture distributed flow section 18.3 (5.58 from the upstream end 
the flume and the surface runoff from the flume are shown Figs. and 
11, respectively. Due very few and scattering experimental points difficult 
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FIG. 10.—Hydraulic Properties Soil Layers Simulated Verification Conjunctive 
Model 
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make comparison for the soil moisture profiles. However, the agreement 
between the observed and computed surface hydrographs reasonably good. 


physically based mathematical model developed simulate conjunctive 
surface-subsurface flow system. Verification the model given comparing 
the numerical results with some existing analytical solutions and limited available 
experimental results. The flexibility and the generality the model allows its 
potential applications quantitative studies variety engineering problems 
including the following which either were already investigated (2), are 
progress planning stages. 


The amount surface runoff, interflow, and ground-water flow resulting 
from rainstorm can determined. 

The validity simple infiltration formulas widely used practicing 
engineers can checked under variety naturally occurring flow situations. 

Hydraulics some surface irrigation methods can studied. 

Ground-water recharge from surface supplies can determined. 

The variations surface flow affected the underlying porous 
medium characteristics can studied. 

The influence the unsaturated zone free surface ground-water flow 
can evaluated. 

Steady and unsteady seepage from wide channels can determined. 

The validity Dupuit-Forcheimer approximation variety flow 
conditions can examined. 
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The following symbols are used this paper: 


S, 


coefficient; 

time rate infiltration; 

Darcy-Weisbach resistance coefficient; 
gravitational acceleration; 

piezometric head; 

index for surface flow sections, also, index for subsurface flow element; 
index for subsurface flow element; 

relative permeability; 

saturated permeability direction; 

saturated permeability direction; 

ground surface hydraulic roughness measure; 
time stage; 

porosity; 

capillary pressure head; 

Reynolds number viscosity; 
time rate lateral inflow; 

degree saturation; 

friction slope; 

bottom slope; 

time; 

volume subsurface flow element; 
cross-sectional average velocity; 

horizontal distance; 

surface flow depth; 

vertical distance; 

time increment; 

surface flow reach length, also, length subsurface flow domain 
element; 

height subsurface flow domain element; and 
weighting factor. 
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TECHNICAL NOTES 


SCALES FOR CIRCULAR JETS 
Nallamuthu ASCE and Thimmaiah 


INTRODUCTION 


Circular jets injected into cross-flows have been studied because their 
application the disposal effluents into rivers and the atmosphere, and 
also because the resulting flow rather complex form turbulent shear 
flow. Consider circular jet diameter, issuing with uniform velocity 
into free-stream large extent, which has uniform velocity U,, 
shown Fig. Experimental observations show that close the jet nozzle, 
due the pressure field imposed the free-stream cross-flow, the jet 
gets deflected. Since the outer layers the jet lose part their momentum 
the entrained fluid due the turbulent mixing, they are more easily deflected 
than the central region, and the jet acquires characteristic kidney shape within 
short distance from the jet nozzle (1). This deformed cross section the 
jet houses its sides pair attached vortices (3,5,8,10). Entrainment 
the free-stream fluid (and thus, momentum) causes further deflection the 
jet, and increases the circulation its cross section, which reflected the 
growth the counter-rotating vortex pair. After some more distance along 
the jet, the jet moving only small differential angle from the free-stream, 
and the cross section the jet almost completely occupied the vortex 
pair. addition this bound vortex system, the cross-flow creates another 
vortex system which shed from the deformed jet (7,8). 

Considering the deflected jet the fully-developed flow region (i.e., after 
the end the potential core), the plane (Fig. 1), the distribution 
the velocity component, (normal the plane) has been found 
similar the direction Keffer and Baines (6) for and 
which the ratio the jet free-stream velocity. The velocity scale was 
taken (u,, U,) which u,, the velocity the axis the direction, 
and the length scale was taken which (1/2) (u,, 
This length scale was found increase linearly with measured from the 
virtual origin, located about 1.3 behind the nozzle, with the coefficient 
proportionality equal 0.215. The velocity scale (u,, U,) was found 
decay much faster than simple circular jet. Based the results 
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FREE STREAM 


FIG. 1.—Definition Sketch: Circular Jet 


experimental study circular jets cross-flows for from 2.7 about 23, 
this note presents practically useful correlations for the velocity and length 
scales. 


AND 


The experiments were performed rectangular open channel with the jets 
issuing from plate fixed parallel and slightly above the bed. For the first 
three experiments with 2.73, 4.52, and 7.05, the diameter the jet was 
0.5 in. (12.7 mm), and the width the flume was in. (22.9 cm), whereas 
for the fourth experiment with 23.4, was equal 1/8 in. (3.2 mm), 
and the width the flume was increased in. (45.7 cm). The depth 
flow was maintained in. (25.4 cm) for all four experiments. 

The axial velocity, the deflected jet was measured with 3-mm external 
diameter Prandtl-type pitot-static tube alined with the The velocity field 
was measured for about seven for the first experiment; for the other 
three experiments, the corresponding values were approx 11, 17, and 
58, wherein measured from the nozzle. More details regarding these 
experiments are available report the writers (12). 

the region developed flow, the plane, the velocity profiles were 
tested for similarity plotting U,)/(u,, U,) versus which 
u,, the maximum value occurring the axis the jet; and 
where 1/2 (u,, U,). was found that these profiles are 
similar for the values tested. was also found that the velocity profiles 
the corresponding length scale (12). These similarity profiles are described 
well the familiar exponential relation (11). 

For these similarity plots, the velocity scale (u,, U,), and the length 
scales are and Considering the maximum velocity, u,,, the results 
the writers and several others, when plotted with(u,, U,)/(U, U,) against 
E/d, showed that the dimensionless excess velocity the deflected jet decays 
more rapidly becomes smaller. Fig. shows one successful correlation 
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developed for the dimensionless excess velocity, which the intercept 


the straight lines the (u,, U,)/(U, U,) versus (double log) 
increases linearly with log-log plot. Thus, using Figs. and (u,, 
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FIG. 4.—Length Scales 


U,) can found for any value within the range studied. has 
also been found Ref. that similar analysis can performed for 
the cross-flow, and empirical correlation for predicting any has been 
developed Ref. 12. 
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Regarding the length scales, shown Fig. both the and scales 
appear increase linearly with (the scales with the minus sign refer 
the negative side). For the three smaller values tested, the scale 
much larger than the scale, whereas for 23.4, the scale only slightly 
larger than the scale. course, for the simple jet (into stagnant ambient) 
for which both the scales will the same. Fig. also shows the 
growth rate the length scale for the simple jet for comparison, which 
slower than the scale and faster than the scale the jets cross-flow 
for the three smaller values. 


Based experimental study circular jets cross-flow for 2.73, 
4.52, 7.05, and 23.4, which the velocity profiles with respect the cross-flow 
have been found similar the and directions the fully-developed 
flow region, the length scales have been found grow linearly. empirical 
correlation has also been developed for the dimensionless velocity scale. 
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TECHNICAL NOTES 


SAND-WATER SLURRY EXPERIMENT 


Experimentally-determined results the fiow sand-water slurry 
8-in. diam (200-mm) polished steel pipe are reported herein. 


Experimental results two tests the 8-in. (200-mm) slurry test loop 
the Georgia Iron Works (GIW) Hydraulic Laboratory are presented (one 
test with water, Test 102-79; the other test with sand-water slurry, 
Test 97-79). Because the modern instruments used these tests and because 
the data-acquisition system, the data are sufficiently precise warrant physical 
interpretation the small difference revealed between the flow water and 
the flow sand-water slurry. 

The principal reason for accuracy the experiments the 32-channel 
data-acquisition and the data-processing system installed the GIW Hydraulic 
Laboratory. All the measuring instruments transduce the sensed physical 
quantities into ma-20 electrical signals which are then processed through 
analog-to-digital converters for further processing through minicomputer. 
means this system, the output all temperature, pressure, and discharge 
measuring instruments are recorded during the same time period. Mean values 
all measured flow variables are displayed cathode ray tube (CRT) for 
the test engineer. The period during which mean values are calculated the 
computer can selected the test engineer. addition, one column 
the CRT display devoted three-digit coded display indicating the change 
value from the previous mean that the test engineer can avoid transients 
the recorded data. When the test engineer satisfied that the flow tolerably 
steady, run taken with the digital data being recorded. 

method insuring accuracy the use redundant instrumentation, 
demonstrated Table Values the discharge, the column labeled 
CH3 were determined means 8-in. Fischer and Porter magnetic flowmeter. 

Emeritus, Georgia Inst. Tech.; and Consultant, Georgia Iron Works Co., 
P.O. Box 626, Grovetown, Ga. 30813. 

Pres. Engrg., Research Development, Georgia Iron Works Co., P.O. Box 
626, Grovetown, Ga. 30813. 
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The values listed the CH6 column were determined means in-situ 
calibrated bend meter. Values the friction-pressure gradient were determined 
from two completely independent sensing and measuring systems for determining 


TABLE 1.—Measured Values and 


i,,, feet water 
gallons per minute per foot 
(2) 


(a) Test 102-79 


unitless but not dimensionless. 


the pressure drop the 50.0-ft (15.24-m) long horizontal, 7.98-in. (203-mm) 
diam, steel pipe test section. 

Sand.—The results sieve analysis the foundry sand, 2.65, are 
listed Table 


Run 
(1) (6) 

280 235 0.000 0.0011 1.00 

482 451 0.002 0.0034 1.00 

715 694 0.006 0.0072 1.00 

877 868 0.010 0.0106 1.00 

1,040 1,013 0.015 0.0150 1.00 

1,219 1,195 0.020 0.0206 1.00 

1,552 1,535 0.030 0.0306 1.00 

1,781 1,755 0.039 0.0397 1.00 

1,991 1,974 0.048 0.0489 1.00 

2,237 2,226 0.060 0.0606 1.00 

2,483 2,478 0.073 0.0735 1.00 

2,713 2,709 0.086 0.0869 1.00 

2,949 2,941 0.101 0.1018 1.00 
3,173 3,164 0.116 0.1163 1.00 

3,406 3,400 0.132 0.1319 1.00 

(b) Test 97-79 
436 470 0.003 0.0042 1.00 
631 601 0.007 0.0083 1.00 
843 819 0.011 0.0121 1.00 
1,000 977 0.019 0.0198 1.00 
1,253 1,197 0.026 0.0269 1.02 
1,413 1,365 0.030 0.0305 1.04 
1,598 1,542 0.036 0.0366 1.04 
1,788 1,726 0.042 0.0425 1.05 
2,061 1,930 0.050 0.0507 1.05 
2,197 2,101 0.057 0.0578 1.05 
2,468 2,384 0.071 0.0714 1.06 
2,746 2,604 0.082 0.0831 1.06 
2,877 2,770 0.092 0.0927 1.06 
3,076 2,967 0.104 0.1052 1.06 
3,470 3,355 0.131 0.1319 1.06 
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The sieve-size data were plotted log probability paper, and straight line 
passed through the 15.9% and 84.1% points from which (geometric mean 
diameter) 270 (geometric standard deviation) 1.6 were determined. 

Test Section.—The 50.0-ft (15.24-m) long horizontal test section 8-in. 
(200-mm) diam steel pipe fitted with two 1/8-in. (3-mm) diam piezometers 
the cross sections the ends the straight test section. The piezometers 
are positioned 45° either side top dead center. The test section preceded 
60.9-ft (18.6-m) length straight horizontal pipe. Each piezometer 
connected means plastic tubing enclosed sediment trap which, 
turn, connected one side Rosemount pressure transducer, thereby 
forming two independent pressure-measuring systems; see and 13, 
Table 

Pressure-transducer calibration checked intervals means open- 
column manometer and water-mercury differential manometer. All pressure 
transducers are calibrated feet water having the specific weight 62.4 
(9802 

Discharge.—The discharge, metered means 8-in. Fischer and 
Porter magnetic flowmeter and bend meter. Calibration both meters 


TABLE 2.—Size Determination 


Sieve number 


Amount finer than, percentage 


(2) 


checked intervals means water flow through standard square-edged 
orifice and volumetric tank. The bend meter utilized backup meter 
because the pressure differences across the bend become very small low 
flows, and because the bend meter questionable instrument when slurry 
containing settling solid forced through the bend. 

Temperature.—Slurry temperature measured means temperature 
transducer cemented the exterior the steel pipe downstream from the 
pump and upstream from the test section. 

Concentration.—Solids concentration determined means pressure 
measurements the upward and downward flowing legs inverted tube. 
The 20.7-ft (6.3-m) high inverted tube instrumented with piezometers with 
8.00-ft (2.44-m) axial spacing both legs. Since measurements are made 
vertical sections pipe, all the solid particles are suspended with the 
result that concentration delivered concentration transport concentration. 
The accuracy the instrument has been checked with the flow water, 
means the Marcy meter with the flow nonsettling slurries and total 
diversion sampler. Based upon the basic theory and the numerous confirmations, 
the values slurry specific gravity, listed the sixth column Table 
are believed accurate the listed significant figures. 


(1) 

96.5 

83.1 

51.4 

25.8 

100 10.3 

140 4.5 
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The concept pseudofluid employed the following. applied 
water-sand slurry, the fluid density, p,,, the mixture density based upon 
the delivered concentration and the fluid viscosity, that the fluid carrier, 
water. 

The flow characteristic greatest interest engineers resistance flow. 
For flow fluid (or pseudofluid), resistance flow described the 
Darcy-Weisbach boundary drag coefficient, For the flow Newtonian fluids 
and the relative roughness, expressed the Colebrook function 


9.35 


which signifies The value determined from experimen- 
tally-determined values the friction-pressure gradient, i,,, and the Darcy-Weis- 
bach equation, follows 


which pipe diameter; pressure gradient horizontal pipe; 


mean velocity; and p,, density the mixture. 

The friction-pressure gradient given Table terms i,,, which has 
the units foot water per foot pipe meter water per meter 
pipe. Conversions are follows 


for pressure gradient having the units pascals per meter. 

the GIW Hydraulic Laboratory, characteristics the flow settling 
slurries are compared the flow characteristics water through the 50.0-ft 
(15.24-m) long, 7.98-in. (203-mm) diam test section. The water test, Test 102-79, 
was performed two days after the slurry test, Test 97-79. The value the 
relative roughness, the test section evaluated from the water test 
finding the value obtain close agreement between the value 
determined from experimental measurements (see Col. Table 3), and 
the value calculated from the Colebrook function (see Col. Table 
3). this procedure, the value was found represent the 
relative roughness the steel pipe the test section the time Tests 
102-79 and 97-79. 


(2) 
(3a) 
for pressure gradient having the units pounds per square foot per foot 
and 
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Test 97-79 enough foundry sand was placed the slurry loop produce 
slurry specific gravity, s,, 1.06, during Runs 11-15, (see Table 1), when 
all the sand was carried suspended load and none bed load. With 
lesser values and more the sand was transported slower-moving 


TABLE 3.—Boundary Drag Coefficients 


Temperature, 
degrees foot per 


(5) (6) 


4 


Note: fps 0.305 m/s; Col. from Eq. Col. from Eq. 


bed load, evidenced the lower values s,,. During Runs and (see 
Table 1), the value s,, 1.0 means that sand was being transported 
either bed load suspended load, i.e., all the sand was deposited 
stationary bed throughout the system. Thus, Test 97-79 encompassed the 


(a) Test 102-79 
84.2 1.79 1.4 0.0172 
84.1 3.09 2.4 0.0157 
84.1 4.59 0.0147 3.5 0.0147 
84.1 5.62 0.0144 4.3 0.0143 
84.3 6.67 5.1 0.0139 
84.3 7.82 0.0144 6.0 0.0136 
84.6 9.95 0.0132 7.6 0.0132 
84.6 11.42 0.0130 8.8 0.0130 
84.9 12.77 0.0128 9.8 0.0128 
85.5 14.35 0.0126 11.0 0.0127 
85.6 15.93 0.0124 12.3 0.0125 
85.9 17.41 0.0123 13.5 0.0124 
86.2 18.92 0.0122 14.7 0.0123 
86.7 20.35 0.0121 15.9 0.0123 
86.9 21.85 0.0119 17.1 0.0122 
Test 97-79 
83.8 2.80 0.0228 2.1 0.0160 
84.0 4.05 0.0216 3.1 0.0150 
85.2 5.41 0.0175 4.2 0.0143 
85.6 6.41 5.0 0.0140 
86.1 8.03 0.0174 6.3 0.0136 
86.9 9.07 0.0153 0.0133 
87.3 10.25 8.1 0.0131 
88.2 11.47 0.0132 9.2 0.0129 
88.9 13.22 0.0118 10.6 0.0127 
89.5 14.09 0.0118 11.4 0.0126 
90.0 15.83 0.0115 12.9 0.0125 
90.0 17.61 0.0108 14.3 0.0124 
90.1 18.46 0.0110 15.0 0.0123 
91.7 19.73 0.0109 16.3 0.0122 
92.3 22.26 0.0109 18.6 0.0121 
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range flow regimes, i.e., stationary bed, heterogeneous suspension, and 
homogeneous suspension. Inasmuch the volume concentration the sand 
was small, i.e., 4%, calculation the boundary drag coefficient, from 
Eq. using the total cross-sectional area, the pipe for calculating velocity, 
reasonable regardless the flow regime prevailing during any run 
Test 97-79. Values calculated from experimental measurements are listed 

order display the effects the sand upon the flow compared 
the writers’ opinion, the results Test 102-79 shown Fig. are indicative 
the veracity the Colebrook function, (see Eq. 1), and the experiment 
that, except for Runs and experimentally-determined values, are 
within two significant-figure agreement with the Colebrook values Eq. 


TEST SYMBOL 
02 - 79 10} 


HETEROGENEOUS HOMOGENEOUS 
SUSPENSION SUSPENSION 


Vv (FPS) 


The experimentally-determined results Test 97-79 are shown solid circles 
Fig. The three flow regimes are labeled the lower portion Fig. 
with the regimes being delineated the deposit velocity and the limit velocity. 

The deposit velocity 8.5 fps (2.6 m/s) was established observing the 
beginning stationary bed the glass section the pipe located downstream 
from the resistance test section. 

The limit velocity fps (4.6 m/s) was established from the values 
the slurry specific gravity, s,,, listed Table During Runs 11-15, Test 97-79, 
1.06. The constant value s,, when fps (4.6 m/s) indicative 
that almost all the sand was carried suspended load with negligible 
portion none being carried slower-moving bed load. 

The anomalous point, Test 97-79 the stationary-bed regime, probably 
results from dunes migrating slowly past the piezometers the test section. 


DEPOSIT 
VELOCITY 
| Limit 
| VELOCITY 
| 
| | 
(0) | | 
STATIONARY 
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Thomas Hagler, Jr., President GIW, encourages the staff the Hydraulic 
Laboratory pursue and report in-house noncommercial studies like the 
one reported herein. John Maffett was the test engineer during Tests 102-79 
and 97-79. 


EQUATIONS FOR SCOUR 
AROUND PIERS 


Christopher James 


INTRODUCTION 


Breusers, al. (2), after extensive study the experimental data 
scour depth around bridge piers sand beds, arrived the following design 
formula 


which the maximum scour depth measured from the mean bed surface; 
the pier diameter; the water depth; the channel mean velocity; 
the value which general sediment motion begins; and 


for (clear water scour); 1.0, 


for (scour with continuous sediment (2) 


Functions and functions the pier shape, incidence the flow, etc., 
which are described (2). Fig. shows the scour depths measured Chabert 


Aerodynamics Section, Railway Technical Centre, British Rail, London Road, 
Derby, England. 

Note.—Discussion open until September 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication July 14, 1980. 
This paper part the Journal the Hydraulics Division, Proceedings the American 
Society Civil Engineers, Vol. 107, No. HY4, April, 1981. ISSN 0044- 
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and Engeldinger (3) around circular piers large variety flow conditions, 
and compares them with the scour depth predicted Eqs. and For these 
tests, 0.67 d,/b 7.0; 33.3 d,/D 1,346; and 0.28 
and the kinematic viscosity. can seen that the predictions are not 
entirely adequate, d,/b being grossly overpredicted some cases. This note 
suggests possible improvements and 


(1) 
After study the mechanics the vortex the scour hole for 


Baker (1) arrived the following formula for scour depth around cylinder 


- @ tan 
b 1 2 


Line of 


slope 
20 2 Error Bounds 


| (Equations 1 and 2 ~ Breusers et al.) 


FIG. Experimental Results Chabert and Engeldinger (3) with 
Formulas Breusers, al. (2) (see Eqs. and 


and d,/D. Eq. can recast the form Eqs. and for 
writing 


which the value which general sediment motion begins (and 
easily shown function d,/D and G); and another function 
and Eq. then becomes 


2.0 2 
: 4 av 
measure 
ft 
“f7 
1.0 2.0 
. . . . . . . . . . . . 
i N. 


TECHNICAL NOTES 


= a, N. 


The form Eq. suggests that Eqs. and might made agree more 
closely with experimental results they are written 


Line of 


/ unit slope 


2.0 


: (Equations 6 and 7) 


FIG. Experimental Results Chabert and Engeldinger (3) with 
Modified Formulas (Eqs. and 


and and functions shape and incidence the flow. 
analysis the experimental results Chabert and Engeldinger 
(3) suggests the following relationships 
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(measured o/@ 
values) 
1.0 
. 1.0 = 
1.625 
G 3 
0.26 log +0.32 for G<355 10°; 
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Variables g,, and can adequately taken equal f,, and 
The scour depths predicted Eqs. and for the experimental conditions 
Chabert and Engeldinger (3) already described, are compared with experimental 
results Fig. can seen that nearly all experimental results lie within 
20% the predicted values. 

Now let consider the physical meaning Eq. can seen that 
high values and (and thus, the scour depth) not depend 
This expected, since function viscosity and thus, indirectly, 
grain-size Reynolds number. high values the Reynolds number (and 
thus, G), viscous effects become unimportant. Similarly, high values 
the scour depth can seen little influenced variations d,/D. 
This again expected, because when (the water depth) very much 
greater than (the grain size), the ratio unlikely influence the scour 
depth. 

Chabert and Engeldinger presented great many experimental results. 
avoid confusion Figs. and simple, but quite justifiable averaging procedure 


TABLE 1.—Comparison with Field Data Melville 


meters 


Tuakau 
Waitangi 
Wanganui 
Matawhero 


has been adopted. Many investigators have realized that for (or 
the scour depth does not depend This reflected Eq. 
Thus, for the results Chabert and Engeldinger, average has been taken 
the experimental measurements scour depth for fixed values 
the number data points displayed Figs. and 

would thus seem that and predict the scour depth more accurately 
than and This some extent confirmed Table where the 
predictions these formulas are compared with the field data Melville (4). 
can seen that the improved design formulas predict the measured depth 
scour rather better than the original formulas for wide range flow conditions. 


Baker, J., Theoretical Approach the Prediction Scour Around Bridge 
Journal Hydraulics Research, Vol. 18, No. 1980, pp. 1-12. 
Breusers, C., Nicollet, G., and Shen, W., Scour Around Cylindrical 


Eqs. 
and 
(see 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
8.7 1.24 3,833 5.82 1.69 6.9 3.6 2.8 
3.2 3.59 3,500 13.4 2.5 4.5 3.4 3.8 
7.4 2.02 431 5,620 2.4 6.9 2.7 4.0 
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Journal Hydraulics Research, Vol. 15, No. 1977, pp. 211-252. 

Chabert, J., and Engeldinger, P., des Affouillements Autour des Piles Ponts,”’ 
Laboratoire Nationale d’Hydrolique, Chatou, France, 1956. 

Melville, W., Scour Bridge Report Number 117, University 
Auckland, School Engineering, Auckland, New Zealand, 1975. 


The following symbols are used this paper: 


oO 


Subscript 


function Eq. 

pier diameter; 
sediment size; 

water depth; 

scour depth; 

functions Eq. 

acceleration gravity; 
functions Eq. 
functions Eq. 
mean channel velocity; 
water kinematic viscosity; 
water density; and 
sediment density. 


conditions start general sediment motion. 
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Errata 


The following correction should made the discussion Diskin (August, 
1980): 


Page 1404, Eq. 19: Should read AH, BH,_,, 
ABG 


an—m 


The writer appreciates the fact that the tabulated Alberta gravel-bed river 
data have been used the discussers, Blaisdell, Thomas, and Agarwal, 
al. test results obtained from laboratory studies and field observations. 

Both Blaisdell and Thomas apply the Alberta gravel-bed river data evaluate 
the state activity the channel bed based initiation motion criterion. 
Blaisdell adopts the Anderson, al. (20) approach which essentially the 
Shields approach with the critical mobility number taken about 0.044. 
Thomas adopts value the critical Shields mobility number 0.060. Neill 
(27) has suggested that value 0.030 adopted for the design stable 
gravel-bed channels. apparent that there universal acceptance 
the value the critical mobility number for gravelbed channels. 

The writer has adopted the Neill (27) criterion for estimating the bed 
gravel-bed channel was motion condition characterized the 2-yr 
flood flow. The Neill criterion based critical average velocity rather than 
critical boundary shear stress the Shields criterion. The Neill criterion 
strictly applicable for relatively uniform gravel-sized materials for two-dimen- 
sional uniform flow and stated as: 

“February, 1978, Aolad Hossain, Ramachandra Rao, and Jacques Delleur 
(Proc. Paper 13565). 

1979, Dale Bray (Proc. Paper 14810). 


Civ. Engrg., Dept. Civ. Engrg, Univ. New Brunswick, P.O. Box 4400, 
Fredericton, N.B., Canada E3B 5A3. 
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which the critical average velocity for initiation motion the gravel 
bed. When the Neill criterion was applied the Alberta data, 16% the 
gravel-bed river reaches show indicated movement; whereas, 27% show indicated 
movement the Anderson, al. criterion tested Blaisdell. 

The difficulty with most the initiation motion criteria that they are 
generally derived from two-dimensional flow straight laboratory channel. 
The plan form natural river not straight and the cross-sectional shape 
varies along the study reach; therefore, the direct applicability any initiation 
motion criterion natural gravel-bed channel estimate best. 

Blaisdell asked about the effect using data from river reaches for which 
has been estimated that the bed material motion for relatively high 
in-bank flows. The writer recognizes that few cases, the bed material transport 
may appreciable the 2-yr flood flow. Bray (2) reported that hydrophone 
data obtained few the river reaches indicated movement material 
over portion all the river bed flows near the 2-yr flood flow. However, 
the data not indicate the bed was motion there was transport 
fine material over the armor layer the gravel bed. The writer considers 
that, for the generalized analysis, the effects bed material transport not 
significantly influence the results presented the original paper. the other 
hand, the writer recommends that the original equations not applied 
gravel-bed river reaches for conditions which result known major bed material 
transport. summary, the channels were straight rectangular channels 
constant bottom slope, more definitive answer could given the questions 
raised Blaisdell. 

Thomas asks the gravel and boulders are exposed all stages the flow 
for the study reaches. far the writer knows, none the river reaches 
experience major modification due large injection fines reported 
Simons, al. (22). 

Thomas discusses the variation Manning’s with stage. The writer stages 
that the average velocity equations the paper are for relatively high in-bank 
flows. The basic Alberta gravel-bed river data (12) confirm that typically 
higher for low flows. The writer considers that average velocity equations for 
relatively low flows for gravel-bed rivers are difficult develop due the 
problem defining average reach properties and due the fact that the flow 
sometimes becomes highly three-dimensional for low flows result the 
jet and wake type flow. The writer recommends that the original equations 

Agarwal, al. state that Lacey’s regime type equation for average velocity 
(Eq. not expected apply gravel-bed channels. The writer would like 
point out that Lacey did incorporate some data from gravel-bed and boulder-bed 
river reaches the analysis from which Eq. was developed. These data 
are shown Fig. the 1946 paper Lacey (15). Thus, does seem reasonable 
expect that the Lacey equation (Eq. should apply gravel-bed river 
reaches for the case relatively high in-bank flows. 

Agarwal, al. apply the Alberta gravel-bed river data test the dimensionless 
flow relation (Eq. 45) developed Garde and Ranga Raju (24,25). Since 
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and the relationship presented Fig. are constant the gravel range; 
and are essentially constant; and the slope the line Fig. 
equal 1/2, then the Ranga Raju relationship presented Fig. can 
rewritten the more common form as: 


The exponent the same that obtained rewriting the Strickler-Man- 
ning equation (Eq. conjunction with Eq. 1). The form Eq. similar 
that Eq. which was developed from best-fit analysis the Alberta 
gravel-bed river data; however, the exponent Eq. 0.28 which significantly 
different from the value 1/6 Eq. 47. 

Finally, additional information, the writer (26) has shown for specified 


assumptions that the optimum value the Keulegan equation (Eq. 


26. Bray, I., Effective Boundary Roughness for Gravel-Bed 
Canadian Journal Civil Engineering, Vol. No. June, 1980, pp. 392-397. 

27. Neill, R., Initial Movement Coarse Uniform Journal 
Hydraulics Research, Vol. No. 1968, pp. 173-176. 


Errata.—The following corrections should made the original paper: 


Page 1104, Table Col. line 10: Should read instead 
Page 1112, the equation under Table Should read 


Closure Gary Parker,’ ASCE 


The discussions Thorne, Alonso, and Bathurst, and Kondap and Garde 
raise interesting and worthwhile points. Thorne, Alonso, and Bathurst raise 
two issues. First, they suggest that explanation nonuniform lateral 
distributions boundary shear stress rectangular flumes based the concept 
They offer alternative explanation terms interacting bank and bed shear 
layers. Most the difference opinion here seems matter semantics, 
since shear stress and momentum flux are the same thing, regardless whether 

1979, Gary Parker (Proc. Paper 14841). 


Assoc. Prof., St. Anthony Falls Hydr. Lab., Univ. Minnesota, Mississippi River 
3rd Ave., S.E., Minneapolis, Minn. 55414. 
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not they are induced molecular turbulent processes. The writer did 
emphasize lateral turbulent diffusion opposed molecular diffusion, for 
obvious reasons the context gravel streams. However, lateral molecular 
diffusion can invoked the laminar case explain the distribution Fig. 

The second point raised Thorne, Alonso, and Bathurst most important 
one. concerns the difficult problem straight-channel secondary currents. 
These currents are indeed much weaker than the primary flow. However, the 
main motivation for neglecting them was the lack satisfactory way 
include them; the analysis can hardly said complete until this done. 
The discussers are correct emphasizing them. 

Recently Ikeda (34) has conducted experimental study self-formed, 
mobile-bed, straight channels nonsuspendable material, which casts some 
light the problem. certain experiments parallel bed ridges occurred, providing 
evidence paired cellular vortices. Measurements downstream mean flow 
velocity distributions revealed that the isovels oscillated the lateral manner 
consistent with the upwelling slow water and downwelling fast water 
that would produced such cells. Ikeda’s data, however, provide strong 
support for Eqs. 21, 22, and 26, even though these equations are based 


Pair vortices 


\ \ 


Sand ridge 


FIG. 10.—Schematic View Inferred Secondary Vortices and Bank Geometry (Cour- 
tesy Ikeda) 


analysis that neglects secondary current. concludes that the cellular nature 
the secondary currents causes their effect rather local that they 
not affect significantly the gross balance that maintains stable banks and 
mobile bed. 

indication the secondary flow structure near the bank, and its effect 
the toe the bank, given Fig. 10, courtesy Ikeda. Some steepening 
the toe the bank apparent. 

Kondap and Garde have provided set dimensionless hydraulic relations 
for sand-bed streams that are the same format Eqs. 29a, 29b, and 
and possess some interesting similarities. comparison the width relations 
the two sets reveals oft-noted but remarkable feature. Consider gravel 
stream with rather typical pavement and sand-bed stream 


with rather typical 0.35 mm; let and their respective widths. 


G 


\ \ 7 
Deformed bed 
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which the widths are meters; and cubic meters per second. For 
the range 500/m’/s, varies from 0.90-1.01. The implication 
that sand-bed streams, the mechanics which differ considerably from 
streams, nevertheless maintain widths comparable their gravel 
counterparts when dominant discharges are also comparable. mechanistically 
well-founded explanation this would invaluable the river engineer. 


34. Ikeda, S., Straight Channels Non-Cohesive Bed Materials. Part 
Technical Report, Saitama University, Urawa, Japan, 1980. 


Errata 


The following corrections should made the discussion Kundzewicz 
(December, 1980): 


Page 2073, paragraph line Should read prove that the heuristic condition 
(Eq. 19) not generally instead prove, that the heuristic condition 
(25) not generally 

Page 2073, paragraph line Should read constraint (Eq. 19) 
instead ‘‘the constraint (25) 


The clear definition and delineation different forms Log Pearson type 
(LP) distribution Fig. the author would welcome many practicing 
engineers and hydrologists together with the new methodology proposes. 
The writer applied the author’s methodology the ALE data six 
rivers Africa, Asia and United Kingdom which had previously applied 


“December, 1979, Erwin Weinmann and Eric Laurenson (Proc. Paper 15057). 


1980, Donthamsetti Veerabhadra Rao (Proc. Paper 15391). 


Engr.-Hydro., Sir William Halcrow Partners, Consulting Engrs., Burderop Park, 
Swindon, Wiltshire SN4 0QD, England. 
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the United States Water Resources Council (WRC) methodology (2,5) and 
observed very good agreement between the two results. example this 
comparative study, the flood magnitudes and return periods for Ethiopian 
river with ALE data are given Table Although both the author’s 
and the WRCs methodologies estimate the parameters the method 
moments, the former, which fits real data LP, looks much simpler 
the writer than the latter, which fits logarithmic data (Pearson type 
distribution). Also, converting the data into dimensionless variables 
X,/X (p. 869), the author’s methodology reduces one the three 
parameters unity. 


TABLE Magnitudes and Return Periods for Ethiopian River with 
ALE Data 


Flood Flow, cubic meters per second 


“From Ref. 


Benson (6), who was Chairman the Work Group Flow-Frequency 
Methods, describes how the choice flood-frequency 
estimating was made after critical comparison six methods: 


single method testing the computed results against the original data 
was acceptable all those the Work Group, and the Statistical 
Consultants could not offer mathematically rigorous method. appeared, 
consequently, that choice could not made solely statistical 
grounds, choice administrative grounds, for which compelling reasons 
existed, was justified. 


was made clear the Work Group that the choice was not 
and study flood-frequency analysis and improvement revisions 
methods’’ was recommended (6). result this, the original guidelines 
(5) have undergone two revisions 1976 (2) and 1977 (10). surprising 
that neither Bobee’s (1,8) important work with respect the mathematical 
properties nor the investigation Matalas, al. (12) the usefulness 
the coefficient skewness regional basis have been incorporated 
these revisions. Other significant contributions similar the author’s have also 
been made since the last revision the WRC guidelines (7,9,11). All this work, 


Return period, years Author 

(1) (2) (3) 
514.3 
598.3 588.1 
661.3 649.7 
749.0 738.2 
820.0 812.1 
100 895.6 892.3 
200 976.5 980.3 
500 1,092.7 1,108.8 
1,000 1,188.5 1,216.4 
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the writer believes, calls for further review and updating the WRC guidelines 
early date. 


Lastly, subsequent paper the same topic (13), the author presents 
the method mixed moments The writer would like know how 


this method compares with the computationally more difficult maximum likelihood 
approach LP. 


Uniform Technique for Determining Flood Flow Bulletin No. 15, 
Water Resources Council, Washington, D.C., Dec., 1967. 

Benson, Flood-Frequency Estimating Methods for Federal 
Water Resources Research, Vol. No. Oct., 1968, pp. 891-908. 

Bobee, B., ‘The Log Pearson Type Distribution: The T-Year 
Event and Its Asymptotic Standard Error Maximum Likelihood Water 
Resources Research, Vol. 15, No. Feb., 1979, pp. 189-190. 

Bobee, B., and Robitaille R., Use the Pearson Type and Log Pearson 
Type Distributions Water Resources Research, Vol. 13, No. Apr., 
1977, pp. 427-443. 

Condie, R., Log Pearson Type Distribution: The T-Year Event and Its 
Asymptotic Standard Error Maximum Likelihood Water Resources 
Research, Vol. 13, No. Dec., 1977, pp. 987-991. 

“Guidelines for Determining Flood Flow Bulletin No. 17A, Water 
Resources Council, Washington, D.C., June, 1977. 

Kite, W., and Frequency and Risk Analysis Hydrology, 
Water Resources Publications, Fort Collins, Colo., 1977, Chap. 15., pp. 

Matalas, C., Slack, R., and Wallis, R., Skew Search 
Water Resources Research, Vol. 11, No. Dec., 1976, pp. 815-826. 

Rao, V., Pearson Type Distribution: Method Mixed Journal 
the Hydraulics Division, ASCE, Vol. 106, No. HY6, Proc. Paper 15477, June, 
1980, pp. 999-1019. 


Discussion Subhash ASCE, and Nevzat Yildirim* 


appears the writers that the shape the assumed temperature profile 
the authors’ numerical experiments was guided the measured temperature 
profiles their laboratory experiments. Though the assumed profile may 
good approximation for very short sidearm, was the case their 
laboratory experiments, not valid for long sidearm, particularly the 
outflow (bottom) layer where the temperature assumed constant and equal 
the temperature along the bottom. According the assumed temperature 
profile, the temperature the water sinking from the inflow layer the outflow 
the following sections. 


“May, 1980, Terry Sturm and John Kennedy (Proc. Paper 15398). 
Engr., Inst. Hydr. Research, Univ. City, 52242. 
“Grad. Student, Civ. and Environmental Engrg., Univ. City, lowa 52242. 
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the authors’ numerical results shown that the velocity and temperature 
approximately decrease linearly and parabolically, respectively. Based dimen- 
sional reasoning the senior writer predicted the same variations for the velocity 
and temperature (11). According his analysis the velocity and the temperature 
distributions long sidearm are given 

342 


which and constants; K/p,c,; surface heat loss coefficient; 

normalized similarity function; prime denotes the differentiation with respect 

Substitution for from Eq. into Eq. and integration yields 


The equations motion (Eqs. and can now written terms 
the similarity functions and Eliminating the pressure term from 
and cross-differentiation, and substituting and from and 
58, found that 


+ 


deriving Eq. 59, the shear stress Eq. expressed 
which k,, turbulent momentum diffusivity. The energy equation (5) becomes 


Eq. differentiation with respect yields 
59, 60, and may written for the special case the condition 
the bottom. both the and velocity components and are 


zero and also the temperature gradient zero there heat flux 
through the bottom boundary; 59, 60, and become 


(f” A A 


eliminating between Eqs. and and solving for one 
gets 
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(f = m 


One can see from that since and then from 
Eq. that One would expect that the temperature would minimum 
the bottom where not zero. Therefore cannot zero anywhere. 


11. Jain, C., Currents Sidearms Cooling Journal the Hydraulics 
Division, ASCE, Vol. 106, No. EY1, Proc. Paper 15320, Apr., 1980, pp. 9-21. 


CHARACTERISTICS Low 
Discussion Michel 


This paper most important attempt summarize the state the art 
and the needs low flow analysis. far frequency analysis concerned, 
clearly emphasizes the urgency developing alternatives the M-day, 
low flow approach. The method proposed Weisman (36) promising. Actually, 
the practical importance characterizing low flows the duration negative 
daily flow runs has been recognized for many years some institutions. 
France for instance, the national applies the concept 
depth-duration-frequency daily flow series, the depth corresponding then 
the threshold flow (40). This partial duration series approach more convenient 
than the annuai maximum series approach this sense that the latter can 
applied only high threshold flows very periodic streams, both, ensuring 
least one negative run per observed year. the other hand, the annual 
maximum series approach may more practical for the fitting theoretical 
distribution. 

Renewal theory teaches that the survivor function R(x) duration, that 
its probability reach larger than uniquely determined 
its age-specific failure rate and has the form: 


1980, the Task Committee Low-Flow Evaluation and Needs the Committee 
Surface-Water Hydrology the Hydraulics Division (Proc. Paper 15400). 

Engr., Hydr. Lab., Dept. Civ. Engrg., Federal Inst. Tech., CH-1015 
Lausanne, Switzerland. 
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However, the verification variable failure rate, that departure from 
the exponential distribution, may lead erroneous conclusions and the use 
physically poorly founded models. Actually, Eq. supposes homogeneous 
process, whereas the succession positive and negative daily flow runs 
generally strongly seasonal phenomenon. the case small streams, the 
seasonal effect often sufficient explanation for departure from the 
exponential distribution; the writer (41) derived the survivor function applicable 


river Rotenbach 

Plaffeien 
record period 1962 -1975 
Qth 
DMAX =74 days Exp. 


0.60 


X=D/DMAX 


station Wasen 


record period 
Qth= 0.00207 


DMAX days 


0.40 
X=D/DMAX 


FIG. 4.—Empirical and Theoretical Log Survivor Functions Negative Runs Duration 
Duration, DMAX Maximum Duration, Oth Threshold 


Wits, 
0.00 0.20 0.80 1.00 1.20 
Eq.2 
Exp. 
0.20 0.80 1.20 
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this case: 


= 


which the cyclic occurrence rate the downcrossings; the 
cyclic failure rate the negative runs duration; and the fundamental period 


river Loire 

Blois 

record period 1863-1974 
Qth= 

DMAX days 


0.40 0.60 


X=D/DMAX 


0.40 0.60 
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FIG. 5.—Empirical (Left) and Theoretical (Simulated, Right) Log Survivor Functions 
Positive Runs Duration 
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(365 days). Fig. shows two examples the fit Eq. two small streams 
Switzerland the maximum likelihood method. 

However, the frequency analysis the duration negative runs has also 
its limitations: (42) most pertinently points out, single negative 
run of, say, days then considered more severe event than sequence 
of, say, 6-days, 9-days, 5-days and 2-day negative runs separated very short 
positive runs. 

Actually, negative and positive runs should analyzed simultaneously. The 
simulation alternating inhomogeneous renewal processes (41) permits 
exhaustive use the observed sequence and downcrossings. For larger 
streams usually necessary introduce into the model. This 
can easily achieved simulating first homogeneous nonpoissonian renewal 
process and then transforming into inhomogeneous process alternating 
rescaling the time coordinate. 

important notice that with low threshold flows, the positive runs 
duration more strongly influenced the seasonal effect than the negative 
one. Fig. shows example where any usual theoretical distribution would 
fail fit the data closely, whereas the seasonal model (transformed gamma 
alternating process) does quite satisfactorily. 


40. Miquel, J., and Phien Bou Pha, B., des Ressources 
des Etiages Estimation leur Probabilité Proceedings the 
Eighteenth Congress the IAHR, Vol. Subject Cagliari, Sept., 1979, pp. 287-294. 

41. North, M., Simulation des Processus Hydrologiques Intermittents par des Modéles 
Alternés Hydrological Sciences Bulletin, Vol. 25, No. Mar., 1980, 
pp. 5-12. 


Schroder, M., Belke, D., and Brandt, Th., der 
Wasser und Boden, No. 1978, pp. 86-89. 


Discussion James Grant Mac ASCE 


The authors have made interesting analysis the Manning Equation 
roughness coefficient natural channels. The roughness coefficients for each 
subsection the channel cross section are determined from velocity measure- 
ments and composite value computed from the average friction slope. These 
values were then compared with composite roughness coefficients computed 
directly from four formulas. 

The results the study found that the composite roughness coefficients 
computed with the Lotter formula have the best correlation with the overall 
value based the apparent friction slope, the writer would expect. The 


“June, 1980, Asok Motayed and Muthusamy Krishamurthy (Proc. Paper 15449). 


Vice Pres., Flaherty Giavara Assocs., P.C., One Columbus Plaza, New Haven, 
Conn. 06510. 
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Lotter equation weighs the roughness coefficients the overall cross section 
proportion the conveyance each subsection. also compatible with 
the method computing the total channel conveyance the U.S. Army 
Corp Engineers HEC-2 computer program. 

The Horton equation for estimating the composite roughness coefficient 
natural channels would not expected match field data because the equation 
does not represent actual flow conditions. The equation assumes that there 
are equal flow velocities all subsection the channel, which condition 
that seldom occurs natural rivers with irregular cross sections. 

The equation Einstein and Banks also simplified expression that does 
not mathematically model actual field conditions. assumes that the hydraulic 
radius each subsection equal that the whole cross section. While 
this assumption reasonable for very wide channels uniform depth, 
not valid for channels with irregular cross sections. 

interesting note that the slope the energy grade line (friction slope) 
varied from one subsection the next when computed from the gaged velocity 
data. This may due either the computation method, the accuracy flow 
measurements, may reflect the effect unknown secondary currents and 
form losses. 

Research the subject composite roughness coefficients very important 
practicing engineers and appreciated those who analyze and design 
channels. would helpful the authors would identify the ranges channel 
widths, depths, and velocities that were used their study, well determining 
the Horton Einstein equations had improved results the wider more 
uniform channel cross sections. 


The boundary layer development region open channel high aspect 
ratio may analyzed use the equation gradually varied flow 


provided that the velocity and pressure distribution correction coefficients are 
close unity. 


1980, Edward Silberman (Proc. Paper 15516). 


Dept. Mechanical Engrg., Univ. Surrey, Guildford, Surrey GU2 5XH, 
England. 


which the rate degradation mechanical energy given 

pgD 
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Whereas fully developed flow, would evaluated from resistance 
formula such that due Manning, boundary layer development necessary 
evaluate the boundary shear, terms boundary layer thickness. 

the case laminar boundary layer 


oy y=0 


and, the turbulent case 


v 1/4 


Provided that some assumption made concerning the origin the boundary 
layer its displacement thickness then Eq. may solved give 

The merit the author’s approach, which fundamentally similar that 
aforementioned, lies the use nondimensional ratios depth and displacement 
thickness most generalized form. Care would, however, need exercised 
the use this approach. 

the two common examples open channel boundary layer 
development are those due flow from sluice gate and over crest. 

The writer has experience boundary layer calculations the latter case 
and would point out that streamline curvature causes free vortex type velocity 
distribution over the crest. 

displacement thickness defined Eq. then, this case, the result 
negative value for because the velocity the edge boundary layer 
greater than that the free surface. 

The definition critical depth given Eq. open criticism. 
clear exposition critical depth given Jaeger (10) terms discharge, 
energy, and momentum and shows that critical depth unique. 

The effect boundary layer distort the velocity distribution, critical 
depth being modified use the usual correction coefficients. 

There can justification for adding displacement thickness the critical 
depth for potential flow. One might just well add the momentum energy 
thicknesses give three different results, which unacceptable. 

The depth defined the author useful reference datum but, though 
close to, not equal the critical depth. 


10. Jaeger, C., Engineering Fluid Mechanics, Blackie and Son Ltd., London, England, 
1961, pp. 110-112 and 138. 


rate degradation mechanical energy; 
boundary layer thickness; 
dynamic viscosity fluid; 
kinematic viscosity fluid; 
density fluid; and 
boundary shear stress. 
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Discussion Subhash Mehrotra’ 


The paper welcome addition the existing literature scour velocities. 
professedly devoted short duration discharges opposed 


a 
W 
W 


USSR CURVE 
CLAY SOILS CURVE 


SILTY CLAY SOILS CURVE 
SANDY SILT SOILS CURVE 


1.0 


FIG. 4.—Correction Factor for Permissible Velocities Function Flow Depth 
long-term constant The writer feels that the distinction between the 
two types flows, far their eroding potentials are concerned, has not 


“September, 1980, William Hughes (Proc. Paper 15701). 
Specialist, Bechtel Inc., Gaithersburg, Md. 20760. 
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been clarified. Are the variable short duration discharges still stage 
evolution and therefore unsteady? the channel velocities refer steady 
states, however short duration? From description given the paper 
the methods employed estimate these velocities, would appear that the 
answer the latter query the affirmative. this the case, the duration 
flow little importance ergodic sense. Thus, the results obtained 
the paper are expected, viz. that accepted maximum permissible 
velocity values appear appropriate for channels subject variable 
intermittent flow conditions well long-term constant flow 

interesting examine the effect flow depth permissible velocities. 
Ref. presents curve from Russian publication, already quoted the paper 
under discussion, which correction factor applied maximum 
permissible velocities plotted against flow depth. The correction factor 
unity for flow depth about 3.2 ft. Fig. presents comparison between 
the curves for clay, silty clay and sandy silt soils derived from Figs. 
and for and the original curve given Ref. Since the derived 
curves are from small-scale figures, they are necessarily subject errors 
interpolation. The agreement among the various curves quite good. This fact 
especially heartening view the fact that the data used Figs. 
and were obtained from two different sources and variety methods were 
employed estimate the channel velocity. 


This paper rich information the characteristics submerged two- 
dimensional horizontal jets over erodible beds. However, deficient its 
use the powerful technique dimensional analysis. result, some 
the conclusions drawn the paper may erroneous, and many the expressions 
derived are very limited applicability. illustrate these points, the authors’ 
investigation the boundary layer thickness, the sand bed critically 
examined herein. 

For given sediment and fluid, proper functional relationship between 
steady-state value and the problem variables can written nondimensionally 
as: 


Table gives experimental values the nondimensional parameters (the first 
three parameters within the parenthesis Eq. 44) derived from the authors’ 


“November, 1980, Chatterjee and Ghosh (Proc. Paper 15820). 
Specialist, Bechtel Inc., Gaithersburg, Md. 20760. 
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data for runs 1-13. Suppose one wishes examine the dependence 
fixed x/L. The proper approach this examination would 
require that all other nondimensional parameters held constant values. 
One must, therefore, choose runs enclosed within bold lines Table only 
for this examination, which and L/D are held constant values (L/D 
actually changes over small range). clear that set such runs 
the range variation u,/ large enough enable one make any 
general conclusion about the relationship between and U,/V gD. Similar 
remarks would apply with respect L/B,. Indeed, the independence 
with respect claimed the authors (Eq. the more surprising, since 
seems intuitively clear that should increase with based results relating 
wall jets. The writer attempted, with considerable strain his eyes, discover 
discernible pattern the scatter experimental points plotted Fig. 
with view detecting any relationship between and with kept constant. 


TABLE Values Nondimensional Parameters for Runs 1-13 


Run number 


(1) 


Comparing run with run 12, for example, could helpful since was 
nearly the same the two runs and varies, though not appreciably. With 
limited number points work with, did appear that indeed increased 
with for fixed U,. (Compare with x.) noted also that because 
large steepness the variation with over the transition and erodible 
bed portions the curve, the actual variation with may quite high 
without appearing so. Thus, the dropping and Eq. appears 
premature and hasty. Only systematic experimental investigation covering 
wide ranges U,/V and L/B, can definitively resolve this question. 
Assuming now that 8/D indeed found independent U,/V 
and systematic investigation the variation with would 
require that done over range constant values The expressions 
derived the authors (Eqs. 16-21) are valid only for one value L/D. 
(Unfortunately even this value changed random, though small amounts, 
the authors’ experiments.) such, these expressions and their nondimensional 


(2) (3) (4) 
0.76 33.0 2.31 
0.92 33.0 2.27 
1.08 2.22 
0.52 22.0 2.31 
0.67 22.0 2.30 
0.95 22.0 2.17 
1.18 22.0 2.21 
0.31 16.5 2.27 
0.59 16.5 2.26 
0.67 16.5 2.24 
0.42 13.2 2.26 
0.49 13.2 2.24 
0.56 13.2 2.22 
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counterparts (Eqs. 22-27) are very limited applicability. 

conclusion, the writer believes that the authors’ study can meaningfully 
extended light the foregoing discussion and hopes that they will undertake 
such extension. The writer also believes that inclusion data scour 
profiles will greatly enhance the value the paper. 
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